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(54) METHOD AND APPARATUS FOR CONTINUOUS CASTING 



(57) Method and apparatus for continuous casting, 
especially casting of steel that can easily provide high 
quality steel that has no central segregation and central 
porosity. In other words, in the method and apparatus, 
central defects are to be eliminated first by identifying 
the solidifying conditions in the full range from the 
meniscus (the surface position of the upper portion of 
molten metal) to the crater end (a final solidification 
position), based on the type (profile) of continuous cast- 
ing machine, type of steel, cross-sectional shape and 
size of a cast piece and the operating conditions such 
as casting speed, casting temperature and cooling con- 
ditions, with special attention paid to the pressure drop 
of liquid phase induced by the liquid flow between den- 
drites resulting from the solidification contraction in 
casting direction in the sol id -liquid coexisting zone, sec- 
ond by calculating the condition of the formation of the 
above internal defects and their positions, and finally by 
applying an electromagnetic body force (Lorentz force) 
in the casting direction in the vicinity of the region where 
the internal defects are formed. 
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Description 

Technical Field 

5 [0001] This invention is concerned with continuous casting, particularly with a suitable continuous casting method and 
apparatus to produce highly qualified steels without segregation and porosity. 

Background Art 

10 [0002] With regard to the continuous castings of carbon steels, low alloy steels, specialty steels and so on, more than 
twenty years have passed since the present vertical-bending-type continuous casting machines began to operate. And 
it has been said that these technologies became mature. On the other hand, the demand for the quality is increasing its 
severity year after year and the pressure to the cost-down also is increasing simultaneously. Aside from the problems 
such as breakout, etc. that often became a problem in the early period of the operating history, there still remains (1) 

15 central segregation and (2) central microporosity as the major problems concerning the quality. 

[0003] The central segregation is the segregation having V characters that takes place with a periodicity in the middle 
of the thickness in the final solidification zone, and is generally called V segregation. The central microporosity is the 
microscopic void that forms in an interdendritic region also in the middle of the thickness in the final solidification zone. 
Summarizing these defects, they are to be called the central defects(internal defects) thereafter in this specification. 

20 [0004] Next, the effects of the central defects on the quality of the steel products are briefly stated as follows. 

(1) The case of thick plate: 

[0005] Hydrogen coagulates and precipitates into these central defects, and so-called hydrogen induced cracking 
25 results during usage. Also, upon welding, the weld cracking occurs starting these defects. 

(2) The cases of rod and wire: 

[0006] Cracking takes place starting the microporosity during drawing. 

30 

(3) The case of thin plate: 

[0007] Upon pressing or during cold rolling, banded defects form, which result from the irregularity in hardness. This 
irregularity is caused by the coexistence of hard and soft spots due to segregation. The above defects takes place dur- 

35 ing the solidification process of continuous casting and lead to a poor quality product. The segregation formed during 
the solidification process remains in final products and can not be eliminated. Tentatively, there is a method of eliminat- 
ing the macrosegregation by diffusion heat treatment. However, this method is not favorable both economically and 
technically because a long period of heat treatment at a high temperature is required. As for the microporosity, it is pos- 
sible to smash them by hot rolling. But whether or not it can completely eliminate them depends on the quantity of the 

40 porosity. Furthermore, an attention must be paid to the fact that the microporosity accompany segregation in many 
cases. 

[0008] Like this, the central defects is the problem associated with the essence of solidification phenomena, and the 
present situation is that it is very difficult to solve by means of the accumulation of know-how or by means of trials and 
errors based improvement. Although there is some difference in degree, these central defects(internal defects) are 
45 common to all the steel grades of slabs, blooms and billets. They exist from the beginning of the continuous casting his- 
tory: Thus, they are an old but at the same time a new problem. 

[0009] Among the measures that have been curried out until now to improve the internal defects, several important 
technologies will be reviewed in the following. 

so (1) Prevention of the bulging 

[001 0] It has been said that the central segregation is formed in slabs with broad width when the solid part of the solid- 
ifying shell or the cast piece between supporting roll pitches was bulged by internal pressure of the steel melt. Although 
this happens by the flow of high solute concentration liquid within the solid-liquid coexisting zone which is induced by 
55 the deformation of the solidifying shell, the detailed mechanism is not clarified satisfactorily. Therefore, to reduce the 
bulging to as small extent as possible, such measures as shortening the roll pitches or dividing one roll into sub-rolls in 
the longitudinal direction have been employed. Besides, the misalignment of rolls is said to be responsible for an inter- 
dendritic fluid flow, thus causing the segregation. However, the internal defects can not be eliminated even if these 
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mechanical disturbances are removed, considering the fact that the central segregation occurs even in the blooms and 
billets where the bulginess hardly become the problem. 

(2) Strengthening of secondary cooling(please refer to Refs. (1) and (2) at the end of this specification) 

[001 1] This is the method of intensively cooling the vicinity of the final solidification position(the crater end) to com- 
press the solid-liquid coexisting phase by contraction force due to thermal stress so as to compensate the solidification 
contraction, thereby reducing central porosity. It has been reported according to the Refs. (1) and (2) that the improve- 
ment was made to some extent. 
io [001 2] On the other hand, the main stream at present is the method of compressing the solidifying shell to give com- 
pressive deformation to central solid-liquid coexisting phase in the vicinity of final solidification position to control the 
interdendritic fluid flow, thereby reducing the internal defects. This method is divided into soft- reduction and hard-reduc- 
tion depending on the amount of reduction. 

is (3) Soft-reduction method at the last stage of solid"rfication(please refer to Refs. (3) and (4)) 

[0013] With this method to improve the central segregation, the solid-liquid coexisting zone is compressed to com- 
pensate the solidification contraction which takes place continuously with the progress of solidification. With respect to 
the soft-reduction amount, a slope needs to be attached so as to correspond to the continuously arising solidification 

20 contraction as precisely as possible. For example, it is shown in Ref. (3) that the central segregation was improved by 
the real machine test of a carbon steel bloom that used the compressive crown roll with roundness attached. Also, in 
Ref. (4), examples are shown about theoretical estimates of reduction gradient necessary for the case of high carbon 
steel blooms (0.7 to 1 wt% C) with 300x500 mm section. According to the estimates, the gradient of 0.2 to 0.5 mm/m is 
required. However, various problems must be overcome to realize this method on a real machine, which will be stated 

25 below. 

® Usually, the soft-reduction is carried out in the range of a couple of meters in the vicinity of the final solidification 
zone, which becomes approximately 0.3 mm/m in the case of the blooms of the above Ref. (4). 

This means that the inclination of 0.3 mm per 1 m needs to be attached to the solidifying shell. Thus, the reduc- 
30 tion amount must be controlled with great accuracy by means of a multi-rolled reduction apparatus, etc. 

@ There is a difficulty that if the amount of reduction is not enough, the effect can not be expected, and that if the 
amount is excessive on the contrary, the interdendritic liquid flows backward to the upstream resulting in the chan- 
nel segregation (i.e. inverse V segregation). 

® Required amount of reduction differs depending on the operating conditions such as a steel grade, dimensions 
35 i: in cross section, casting speed and cooling condition. Therefore, a great amount of labor and cost is necessary for 
these trials and errors to find an appropriate condition even in the case of a few steel grades. 
® Since the soft reduction method often gives rise to the new problem of internal cracking(Ref. (5)), a considera- 
tion must be taken into to prevent this. 

40 [001 4] Thus, it is not easy to make use of this method to achieve the purpose. 

(4) Continuous forging method (Refs. (7) and (8)) 

[0015] Next, stated is the hard-reduction method in which the solid-liquid coexisting phase in the vicinity of the final 
45 solidification zone is mechanically largely deformed thereby squeezing the high solute concentration liquid to the 
upstream to prevent the central segregation (V segregation). There are two methods in this: One is to use large diam- 
eter rolls (Ref. (6)), and the other is the continuous forging method in which the shell is continuously forged using anvils 
(Refs. (7) and (8)). Because both belong to the same category in their concepts, only the latter is described in the fol- 
lowing. As shown in Figure 42, the vicinity of the final solidification zone is forged by anvils while moving toward casting 
so direction together with the anvils. It has been reported that by repeating this cyclically, the high solute concentration liq- 
uid within the solid-liquid coexisting zone is squeezed into the upstream region with low volume fraction of solid, thus 
enabling it possible to suppress the central segregation. Also, it is said that the internal cracking can be avoided by set- 
ting up an appropriate forging condition. It is possible to control the segregation ratio Ke (= ~C/C 0 , ~C = average solute 
concentration. C 0 = solute content) to be Ke <1 by changing the volume fraction of solid at the time of forging. 
55 [001 6] The most important point of this method is to clarify the flow phenomenon in the solid-liquid coexisting zone at 
the time of forging. However, the authors of this reference have derived the relationship between segregation ratio Ke 
and the volume fraction of solid at the time of squeezing taking into account only the conservation law of solute ele- 
ments. In their model, the liquid flow in the solid-liquid coexisting zone has not been treated explicitly, that is to say, the 
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influences of the flow of solute concentrated liquid in the dendritic scale on the segregation has not been clarified. 
Therefore, while it is controllable as for the average macrosegregation in a macroscopic inspection range of the solid- 
liquid coexisting zone, the information about so-called semi-macrosegregation in much smaller inspection range (den- 
dritic scale) can not be obtained. The semi-macrosegregation remains to some degree in their method. 
5 [001 7] Accordingly, the mechanism of the formation of the semi-macrosegregation belongs to a future subject and the 
flow phenomenon of the ejected liquid phase needs be clarified, in this connection, there is a possibility that the V seg- 
regation has already been formed when forged: in this case, the questions are raised on how the ejected liquid behaves, 
on if it remains as the semi-macrosegregation, etc. 

[0018] In the above references, the blooms having nearly square cross section have been studied where the shape 
10 of solid-liquid coexisting zone can be approximated by a cylindrical form, and so when the solidifying shell is com- 
pressed in an iso-concentric fashion, the flow pattern will become comparatively simple. But in the case of slabs having 
broaden width, it is questionable whether or not the flow becomes a simple upstream pattern. In any case, it is not easy 
at all to predict the flow pattern of the solute concentrated liquid and to evaluate its influences when the solid-liquid 
coexisting zone is mechanically deformed to a large extent. 

15 

(5) Electromagnetic stirring (Refs. (9) and (10)) 

[001 9] This is the method of stirring the solid-liquid coexisting zone by an electromagnetic force in the vicinity of the 
final solidification zone to disperse the central segregation (Ref. (9)). For example, there is a method of spiral-stirring 
20 within the cross section of a solidifying shell. Another method is that the electromagnetic force is applied within the sec- 
ondary cooling zone or within the mold with the aim of transiting columnar structure to equiaxed structure. The latter 
method is based on the prerequisite that the equiaxed structure is preferred to the columnar structure as for the central 
segregation, but the theoretical basis is poor. These methods are not an essential solution and are not the mainstream 
at present 

25 

(6) The combination of the above methods (1) to (5) 

[0020] The bulging prevention measure has been esteemed consistently until the present as a fundamental technol- 
ogy and the following combinations are carried out based on this. 

30 [0021 ] For example, it has been reported for carbon steel slabs of 0.08 to 0. 1 8 wt% C in Ref. (10) that with the com- 
bination of short roll pitch with sub-segmented rolls (bulging prevention), taper-alignment method (gradually narrowing 
the gaps between the rolls in the downstream direction to correspond to the contraction of the cast piece due to solidi- 
fication contraction and temperature drop and strengthened cooling + electromagnetic stirring in the secondary cooling 
zone, the central segregation was improved compared to the case with no measures taken. 

35 [0022] It is also stated in Ref. (11) that the central porosity reduces when the equiaxed structures are developed by 
simultaneously adopting low casting temperature and electromagnetic stirring for the carbon steel blooms and round 
billets in which the equiaxed structures are difficult to develop. Furthermore, it is reported that the central segregation 
and porosity can be reduced by adjusting the reduction amount in the final solidification zone and by developing the 
equiaxed grains via electromagnetic stirring within the mold. 

40 

(7) Cast Rolling method in a thin slab casting 

[0023] So-called mini-mill, that concisely sums up a steel making plant, has become increasingly popular because of 
such advantages as recycling of raw materials, the energy saving, a low plant construction cost and the gentleness to 
45 the earth environment in comparison with a large scaled conventional blast furnace. 

[0024] With the mini-mill, thin slabs with the thickness of 50 mm or 60 mm (so-called near-net-shape-castings made 
as close to the size of the final products as possible) are cast, instead of large sectioned conventional castings with the 
dimensions of 200 mm or 300 mm. 

[0025] Here, The Cast Rolling method (Ref. (12)) will be described as an example. This method is characterized by 
so gradually compressing and thinning solidifying shell (reduction ratio of 10 to 30%) by rolls in the range including the 
solid-liquid coexisting and liquid phases. This method is supposed to be born from an idea that the solidifying shell can 
be reduced during solidification considering that there is a limit to make thin at inlet nozzle position, by which the follow- 
ing effects are reported. 

55 (J). Because dendrites are mechanically destroyed, equiaxed fine grains are formed. 
©. As a result, the macrosegregation is fairly decreased. 

[0026] However, the flow behavior of high solute concentration liquid induced during heavy deformation of the solid- 
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liquid coexisting zone is unpredictable, and therefore it is not easy to control so as to avoid detrimental defects such as 
inverse V segregation, etc. 

[0027] So far, key technologies to improve the internal quality of continuous castings of steels were reviewed from a 
vast amount of literature. Historically speaking, they trace back to the taper-alignment method for the control of the bulg- 
ing that causes segregation, progress into the shortened roll pitch/divided roll method, strengthened secondary cooling, 
electromagnetic stirring and presently become soft/hard-reduction methods or the combination of the electromagnetic 
stirring and the soft- reduction. Although the technology has been improving meanwhile, it has not yet reached to the 
essential solution of the problem. 

Disclosure of Invention 

[0028] All of the aforementioned technologies are trials and errors based measures based on the empirical and qual- 
itative insights into the solidification phenomena. Therefore, it is inevitable to take a vast amount of time and labor to 
obtain appropriate conditions every time when the steel grade, the shape and dimensions of the cross section, the 
machine profile and the operating conditions (casting speed, temperature, cooling method, etc.) are changed. Despite 
that, many cases have been seen that the optimal conditions can not necessarily be obtained. In conclusion, although 
the individual technology has succeeded in reducing the segregation temporarily to some extent, there was the incon- 
venience that the essential solution of the problem can not be obtained, because the solidification behavior is not 
grasped based on solidification theory, or more precisely speaking, because the mechanism of the formation of the cen- 
tral defects is not satisfactorily clarified. 

[0029] The purpose of this invention is to solve the aforementioned inconveniences in the conventional technologies, 
and to provide with the method and apparatus especially in the continuous casting of steels that always can get the 
good steel with no central segregation and porosity easily, even if the steel grade, the shape and dimensions of cross 
section, the machine profile, the operating conditions (the casting speed, temperature, cooling method, etc.) are 
changed or furthermore even if the casting speed is increased to raise productivity. 

[0030] Thereupon, this invention tries to eliminate the above-mentioned internal defects by exerting an electromag- 
netic body force (Lorentz force or simply termed electromagnetic force) toward the casting direction in the solid-liquid 
coexisting zone which is prolonged along the casting direction at the central part of cast piece. The aim is to complete 
the feeding of interdendritic liquid toward the casting direction. More specifically, this invention investigates the solidifi- 
cation mode of whole range from meniscus (top surface position of the melt) to the crater end (final solidification posi- 
tion) when the type of the machine, the steel grade, the shape and dimensions of cross section and the operating 
conditions (casting speed, casting temperature, cooling condition, etc.) are given, with the special attention paid to the 
pressure drop of the liquid phase due to the interdendritic liquid flow (Darcy flow) which is caused by the solidification 
contraction in the solid-liquid coexisting zone. This invention possesses the calculation means to figure out the condition 
for the formation of the internal defects, the position of the formation and the electromagnetic body force required to 
suppress the internal defects. And it is equipped with the electromagnetic booster (Exerting means for electromagnetic 
body force) to exert the above-mentioned electromagnetic body force toward the casting direction in the vicinity of the 
position where the internal defects are formed. Thus, this invention is comprised by the above-mentioned constitution 
thereby trying to achieve the aforementioned purposes. 

Brief Description of Drawings 

[0031] 

Figure 1 is the schematic diagram composing the Electromagnetic Continuous Casting system by this invention. 
Figure 2 is the schematic diagrams for describing the details of the electromagnetic booster of Figure 1. 
Figure 3 is the diagrams for explaining solute redistribution of alloy elements, (a) shows the equilibrium phase dia- 
gram for Fe and a certain alloy element, (b) shows the solute distribution for the case of an equilibrium solidification 
type alloy and (c) shows that for the case of a non-equilibrium solidification type alloy. 
Figure 4 shows local linearization model of a nonlinear binary phase diagram. 
Figure 5 is the schematic diagram showing a dendritic solidification model. 

Figure 6 shows the formation site of microporosity and the space size of interdendritic liquid, (a) the formation site 
of the porosity, (b) and (c) the model for calculating the space size of the liquid. 

Figure 7 shows the schematic diagram of the volume element used for the numerical analysis. V L is flow velocity 
vector of the interdendritic liquid and V s is deformation rate vector of the dendrite crystal. 

Figure 8 is the diagram for explaining discretization (quoted from p. 97 of Ref. (20)). Shaded area shows the control 
volume element and the points denoted by the circles are called grid points. Fe, Fw, Fn, Fs at the control volume 
faces e, w, n, s show the incoming and outgoing of a physical quantity <(>. 



5 



EP 0909598A1_L> 



EP 0 909 598 A1 



Figure 9 (a) shows the coordinates system used for the numerical analysis and (b) the topology of the discretiza- 
tion. The meanings of the symbols in (b) are similar as those of Figure 8. 

Figure 10(a) is the outline of the main program that shows the flow of the calculation in the numerical analysis. 
Figure 10(b) is the outline of the flowchart that shows the fluid flow analysis by the momentum equation in the 
s numerical analysis. 

Figure 1 1 shows the numerical results of a large stee! ingot (1 m diam. x 3m height) chosen as an example to show 
the validity of the numerical analysis, is the casting design, (b) shows an example of the contours of temperature 
during solidification (after 11.5 hours). Symbol S denotes solid, symbol M solid-liquid coexisting zone (mushy 
zone), symbol C shrinkage cavity and the broken lines show the boundaries of these phases, (c) is the contour of 
10 the volume fraction of solid at the same time as (b). Numbers in the diagram denote the volume fraction of solid (0 
to 1). (d) shows the liquid flow pattern after 4.28 hours when the whole region became mushy. The portion with high 
density of streamlines means where the flow velocity is high. It is about 3.5 mm/s at the center where the velocity 
is high, (e) shows the interdendritic liquid flow pattern after 1 1 .5 hours. The velocity is about 0. 1 mm/s at the center. 
(0^ and (g) show the distributions of the macrosegregation after solidified. The segregation is expressed by 
15 (C - C 0 /C 0 ) x 100% (C= calculated concentration, C 0 = alloy content), (f) is the case of carbon: the positive seg- 

regation is more than 30% at the central portion and about 5% at the upper O.D. where the A segregation usually 
takes place. The negative segregation is biggest in the lower part near the center (-21%) and becomes lesser near 
middle O.D. (-10%) (g) is the case of phosphorus which shows the same trend as carbon, although both positive 
and negative segregation emerge more prominently. 
20 Figure 1 2 shows the interdendritic liquid flow pattern induced by solidification contraction that was confirmed by the 

numerical analysis (a). The schematic diagram of V defects is shown in (b). Flow velocity is high at the central por- 
tion as denoted by high density of streamlines in (a). The flow to the lateral direction is extremely smaller compared 
to that to the casting direction, (b) shows the V defects that in dendritic scale possess locally prominent positive 
segregation (+) and simultaneously accompany microporosity along the V pattern. The arrows indicate the inter- 
ns dendritic liquid flow where the liquid of the surroundings flows in along the V defects. 

Figure 13 is a schematic diagram of a typical conventional vertical billet caster. Symbol L denotes liquid zone, M 

solid-liquid coexisting zone, S solid zone. 

Figure 14 shows the linearized data of Fe-C phase diagram. 

Figure 15 shows the relationships between the temperature and volume fraction of solid for the steels used by the 
30 numerical analysis. They were obtained by using the nonlinear multi-alloy model, (a) is in the case of 1 C-1 Cr bear- 
ing steel and (b) is in the case of 0.55% carbon steel. 

Figure 16 shows the effects of oxygen content on the equilibrium CO gas pressures in interdendritic liquid phase 
with no CO gas bubbles (refer to eqs. (49) to (58)). (a) is in the case of 1 C-1 Cr bearing steel and (b) is in the case 
of 0.55% carbon steel. 

35 Figure 17 shows the distributions of temperature T and volume fraction of solid gS at the center element (a) and 
the distributions of solidifying shell thickness (b) at the steady state in the best mode No.1 for carrying out the inven- 
tion (vertical continuous casting). The distributions (1 ) show the case that only the temperature was calculated and 
the distributions (2) the case that the thermal conductivity of the liquid was multiplied 5 times considering the Darcy 
flow. 

40 Figure 18 shows the results of the computation No.1 at the steady state of the best mode No.1 for carrying out the 
invention, (a) is the distributions of the temperature T, the volume fraction of solid gS, the liquid pressure P and 
Darcy flow velocity at the center element, (b) is the distributions of the surface heat transfer coefficient H and the 
solidifying shell thickness, (c) is the distributions of the permeability K and the body force X (gravitational / Lorentz 
force) at the center element, (d) is the distribution of surface temperature T s . 

45 Figure 19 shows the results of the computation No. 2 of the best mode No.1 for carrying out the invention. 

Figure 20 shows the phase distribution in the computation No. 1 of the best mode No.1 for carrying out the inven- 
tion. L denotes the bulk liquid region, M the solid-liquid coexisting zone and S the solid zone. The region with more 
than 1% of volume fraction of solid is regarded as M. 

Figure 21 shows the interdendritic liquid flow in the vicinity of the crater end in the computation No. 1 of the best 

so mode No.1 for carrying out the invention. 

Figure 22 shows the dendrite arm spacings in the best mode No.1 for carrying out the invention, (a) is of the case 
using theoretical equations (28) and (29), and (b) is of the case using the empirical equation (31). Since the forms 
of the equations (31) and (71) are different upon the calculation of (b), A and n are determined as A=7.28 and 
n=0.39 respectively so that d becomes d=35 ujti at the surface element where accelerated solidification phenom- 

55 enon is absent. 

Figure 23 is the schematic diagram of the electromagnetic booster in the best mode No.1 for carrying out the inven- 
tion, (a) shows the outlook and (b) the horizontal cross section. The electromagnetic body force (Lorentz force) is 
applied downward in vertical direction. 
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Figure 24 shows the effect of the electromagnetic body force (Lorentz force) in the computation No. 3 of the best 
mode No.1 for carrying out the invention. 

Figure 25 shows the specific heat C (cal/g°C) and the thermal conductivity X (cal/cms°C) of 0.55% carbon steel. 
Figure 26 shows the schematic diagram of the typical vertical bending continuous caster used for the best mode 
No.2 for carrying out the invention. In the diagram, the supporting rolls and water-spray cooling unit are not shown 
except for bending and unbending rolls. 

Figure 27 shows the results of the computation No.1 at the steady state of the best mode No.2 for carrying out the 
invention, (a) is the distributions of the temperature T, the volume fraction of solid gS, the liquid pressure P and 
Darcy flow velocity at the center element, (b) is the distributions of the surface heat transfer coefficient H and the 
solidifying shell thickness, (c) is the distributions of the permeability K and the body force X (gravitational / Lorentz 
force) at the center element, (d) is the distribution of surface temperature T s . 

Figure 28 shows the results of the computation No.2 of the best mode No.2 for carrying out the invention. 

Figure 29 shows the effect of the electromagnetic body force (Lorentz force) in the computation No. 3 of the best 

mode No.2 for carrying out the invention. 

Figure 30 shows the results of the computation No.3 of the best mode No.2 for carrying out the invention, (a) shows 
the solidification profile and (b) Darcy flow pattern in the neighborhood of the crater end. L denotes the bulk liquid 
region, M the solid-liquid coexisting zone and S the solid zone. The distance from the meniscus is that along the 
central axis of the slab. The slab is curved in fact, but is shown in a prolonged rectangular form for the convenience 
of display. The unbending rolls are denoted by circles. 

Figure 31 shows the effect of the electromagnetic body force (Lorentz force) in the computation No. 4 of the best 
mode No.2 for carrying out the invention. 

Figure 32 shows the electric conductivity a (1/nm) of carbon steel (The Iron and Steel Institute of Japan: Iron and 
Steel Handbook 3rd edition, p. 311). 

Figure 33 shows the results of the computation No. 1 of the best mode No.3 for carrying out the invention. 
Figure 34 shows the results of the computation No. 2 of the best mode No.3 for carrying out the invention. 
Figure 35 shows the effect of the electromagnetic body force (Lorentz force) in the computation No. 3 of the best 
mode No.3 for carrying out the invention. 

Figure 36 shows the results calculated to preliminarily examine the effects of soft-reduction in the best mode No.3 
for carrying out the invention, (a) shows the distribution of the reduction required to compensate the net solidifica- 
tion contraction (computed toward downstream from the reference position of 25m from the meniscus where the 
volume fraction of solid at the center is 0.1). (b) shows linear reduction gradients in the neighborhood of the crater 
end. (c) shows the computational results showing the degree of relaxation of the liquid pressure drop for the given 
reduction gradients. 

Figure 37 shows the combined effects of the electromagnetic body force (Lorentz force) and the soft-reduction in 
the computation No.4 of the best mode No.3 for carrying out the invention. 

Figure 38 shows the TTT diagram of 0.55% carbon steel where the symbol A denotes austenite and P pearlite. The 
solid lines indicate experimental data (from Ref. (30)) and the broken lines the calculated values using, eqs. (34) 
and (76) in this description (Ret (21 )). The start and end of the transformation is defined as the volume fraction of 
pearlite becomes g P =0.0l and g P =0.99 respectively. 

Figure 39 is the schematic diagram showing the attractive forces generated between the coils in the case that 
superconductive coils are used as an apparatus to generate DC magnetic field, (a) shows cylindrical coordinates 
system (r, 6, z) . (b) shows the results calculated for the cases that the magnetic flux densities Bz at the center z=b/2 
are set 1,2 and 3 (Tesla) (a is fixed at 0.8m). Symbol I denotes the electric current in coil, and the pressure 
(Kgf/cm2) are the values of the attractive forces divided by the cross-section area of coil. 

Figure 40 is a schematic diagram for explaining the relationship between magnetic attraction and reduction gradient 
(mm/m). Most of the deformation is concentrated to the dendritic skeleton of the central part of solid-liquid coexist- 
ing zone where the mechanical strength is extremely low in comparison with solid, and the relationship between the 
magnetic attraction and the reduction gradient (in other words displacement) becomes slightly nonlinear. 
Figure 41 is the diagrams for explaining the staggered grids used for the discretization of momentum equation, (a) 
is the staggered grid in X,(r) direction, (b) is that in X 2 (z) direction and (c) is that in X 3 (Y) direction respectively. 
Figure 42 is the schematic diagram of conventional continuous forging method. Shown is the flow manner in which 
the liquid phase in the solid-liquid coexisting zone is squeezed by anvils to the upstream. 8 indicates reduction 
quantity. 

Figure 43 shows the formation of centerline shrinkage in cast steel (from Ref. (14), p.242). 

Figure 44 is the schematic diagram of the experimental apparatus for pressurized casting. 

Figure 45 shows the measured and calculated temperature histories in the atmospheric casting experiment. 

Figure 46 is the sketched diagram showing as cast macrostructure in the atmosphere. 

Figure 47 is the microstructure of the V pattern in as cast specimen in the atmosphere. 
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Figure 48 shows the variation in Vickers hardness in the vicinity of the V pattern in as cast specimen in the atmos- 
phere (load 1kgf and 10 sec). 

Figure 49 shows the results of the numerical analysis for as cast specimen in the atmosphere, (a) is the volume 
fraction of solid distribution 55 seconds after pouring when internal porosity begins to form, (b) is the volume frac- 
5 tion of porosity (percent) distribution after solidified. 

Figure 50 shows the macrostructure of the pressurized casting at 10atm. 
Figure 51 shows the macrostructure of the pressurized casting at 22atm. 

Figure 52 shows the effects of the pressurized casting predicted by the numerical analysis. The calculated volume 
fractions of internal defects are shown for the cases: (a) atmospheric casting (no pressurization), (b) pressurized 
10 at 1 0atm and (c) pressurized at 20atm. 

Figure 53 shows the effects of the pressurized casting predicted by the numerical analysis for the steel castings of 
Ref. (34). The volume fractions of porosity are shown for the cases: (a) atmospheric casting (no pressurization) and 
(b) pressurized at 4.2atm. 

Figure 54 is the schematic diagram for explaining the mechanism of the formation of internal defects. 
15 Figure 55 is the outline of the bending type bloom caster used for the best mode No.4 for carrying out the invention. 

Rolls other than the unbending rolls are not shown. 

Figure 56 shows the results of the numerical analysis of conventional casting method for the best mode No.4 for 
carrying out the invention. 

Figure 57 shows the effect of electromagnetic body force in the best mode No.4 for carrying out the invention. 
20 Figure 58 is the specific diagram in which the electromagnetic booster by this invention is installed to the continu- 

ous castings having rectangular cross -sections such as bloom and billet, (a) shows cross-section diagram and (b) 
shows the AA side view. The broken lines of (a) denote DC magnetic field and the arrow of (b) denotes casting 
direction. 

Figure 59 is the ground plan of the electromagnetic booster of Figure 58. (a) shows the BB cross-section of Figure 
25 58 and (b) the race-truck-type superconductive coil. 

Figure 60 shows the connection diagrams of DC electrodes, (a) is parallel type, (b) series type and (c) mixture type. 

Figure 61 shows an example of the load distribution when soft-reduction gradient is given to cast piece. 

Figure 62 shows a situation that a gas shield box is attached to prevent oxidization, (a) shows the side view and (b) 

the ground plan. Symbol 108 denote shows plane milling tool. 
30 Figure 63 is the specific diagram of the electromagnetic booster by this invention in which the distance between the 

superconductive coils is narrowed in comparison with that of the booster of Figure 58. (a) shows the cross-section 

and (b) the horse-saddle type superconductive coils. 

Figure 64 is the specific diagram of the electromagnetic booster by this invention applied to the continuous castings 
having wide rectangular cross-sections such as slab. The diagram is that of the side sectional plan. Symbols 129 
35 and 130 denote upper and lower split rolls respectively. 

Figure 65 is the specific diagram (the cross-section) of the electromagnetic booster by this invention applied to twin- 
type continuous casting. Symbol 131 denotes flexible bus bar or cable. 





Description of Symbols 
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[0032] 




1 


Electromagnetic booster 




1a 


High rigid frame 


45 


1b 


Cast piece 




1c 


DC rotating electrode 




1d 


Spring 




1e 


Fixed axle 




1f 


Nonmagnetic roll 


50 


2 


Ladle 




3 


Tundish 




4 


Nozzle 




5 


Water-cooled mold 




6 


Cast piece 


55 


7 


Bending rolls 




8 


Unbending rolls 




9 


Detecting section 




10 


Computer, CPU 
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Final controlling section 
Display unit 



Symbols 13 to 101 missing 
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102 


Electrode 




103 


Flat bus bar 
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104 


L-shape bus bar 




105 


Insulating electrode box 




106 


Spring 




107 


Electrode-fixing frame 




108 


Plane milling toll 


15 


109 


Gas shield box to prevent oxidization 




110 


Electrode box room 




111 


Plane milling tool room 




112 


Gas inlet 




113 


Gas inlet 


20 


114 


Cutting tool 




115 


Discharaina outlet of cut chiDS 




116 


Air gap to release oxidation-preventing gas 




117 


Upper frame 




118 


Lower frame 


25 


119 


Pillar 




120 


Suoerconductive coil 




121 


Goolina chamhpr fnr ^unprmnHnrtix/P mil 




122 


Rigid frame to hold superconductive coil 




123 


Outer cooling chamber 


30 


124 


Upper rolls 




125 


Lower rolls 




126 


Bearing 




127 


Oil-hydraulic cylinder 




128 


Rigid frame both f ixable and movable in the longitudinal direction of cast piece 


35 


129 


Upper split rolls 




130 


Lower split rolls 




131 


Flexible bus bar or cable 



Best Mode of Carrying Out the Invention 

40 

A. Numerical Analysis of Solidification Phenomena 

[0034] In order to know the position of internal defects and their morphologies precisely, the mechanism of the forma- 
tion of the internal defects must be clarified, and further it is inevitable to perform numerical analyses of the solidification 
45 phenomena on the basis of solidification theory. Then first, the theory of the numerical analysis in the computational 
means for this invention will be described in detail. Subsequently, the mechanism of the formation of the internal defects 
will be stated. 

A- 1. Theoretical Equations for the Numerical Analysis of the Solidification Phenomena 

50 

[0035] The equations are described in the following sections that these inventors conceived for the numerical analysis 
of solidification phenomena on the basis of the solidification theory. 

(1) The Energy Equation 



[0036] The energy equation is given by Eq. (1) that describes the energy conservation for a certain volume element 
in the solid-liquid coexisting zone. As shown in Figure 7, the volume element is regarded sufficiently large compared to 
dendrite arm spacing (spacing between the branches of dendrite crystal) and small enough to be able to judge the 
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changes of such physical quantities as temperature T, volume fraction of solid g s , etc. 

— (cpT) + v\c L P p L g L v L +c s p p sgs v s )T)=V-{XVT) + S {1) 



10 

[0037] The details of each symbol are listed in Table 1 at the end of this description. The first term on the left of the 
equation is the change of the heat capacity per unit volume and time, the second term is the divergence (outgoing heat 
quantity per unit volume and time) due to interdendritic liquid flow and deformation of solid, the first term on the right is 
the divergence due to heat conduction and S heat source term. S consists of the sum of the latent heat of fusion, the 
15 effect of solid deformation and the heat of Joule by electric current as shown in Eq. (2). 

S = pL(^ff- + v s -Vg s ) + Qj (2) 
o t 



[0038] Also, the average heat capacity cp in Eq. (1) is given as Eq. (3) using the volume fractions of solid g s and liquid 

25 g L . 

cp = c L pp L g L + Cpp s g s (3) 

30 [0039] Here, introducing the volume fraction of porosity g v , the following relation holds. 

Qs + Ql+Ov^ 1 (4) 

35 [0040] Also, the temperature dependencies of specific heat C, density p and thermal conductivity X are taken into 
account both for solid and liquid. The suffix L denotes liquid and the suffix S solid. Furthermore, Eqs. (1) and (2) can of 
course be applied to the liquid and solid phases and the phases including porosity in addition to the solid-liquid coex- 
isting zone. 

40 (2) The Solute Redistribution Equation 

[0041] Solute atoms dissolve in solid and liquid, and their distributions are determined by the equilibrium phase dia- 
grams and the diffusion rates in each phase. For example, carbon atoms diffuse promptly in solid phase (at high tem- 
perature) as well as in liquid phase. On the other hand, the diffusion rate of silicon atoms in solid is very slow. 

45 Thereupon, it is assumed in this invention that while all alloy elements diffuse completely in interdendritic liquid phase, 
only carbon diffuses but other elements do not regarding the diffusion in solid. In other words, carbon is regarded as an 
equilibrium solidification type element as shown in Figure 3(b) and others non-equilibrium solidification type elements 
as shown in Figure 3 (c). Considering the nonlinearities of liquidus and solidus lines in equilibrium phase diagram as 
shown in Figure 4, the relationship between the solute concentration of solid 0 s * and that of liquid C L at the liquid-solid 

so interface is expressed by Eq. (5) (please refer to Ref. (15) for detailed derivation). 
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Of Of - A n.k(n)Cn + B n,k{n) (5) 
where, 

1 5 (6) 

m n,k(") 



R _ m nMn) L 

m n,k(n) 

15 

m L and m s are the slopes of liquidus and solidus lines respectively, and other symbols are shown in Figure 4 where the 
suffix n denotes alloy element and the suffix k the locally linearized segment number of the liquidus and solidus lines. 
In order to derive the conservation law of solute elements within liquid and solid phases, it is necessary to take into 
account the flow of solute concentrated liquid and also the deformation in the mushy phase. The solute conservation 
20 law including these effects is expressed by the next equation. 

~(QTp) + ^(p L g L ciy L +P S g s Q\-) = V-(l^p L g L vc n L ) (8) 

25 



[0042] The first term on the left of Eq. (8) is the change in the average solute mass, the second term is the divergence 
30 due to the interdendritic liquid flow and the deformation of the mushy zone, and the right side is the diffusion term in the 
liquid. The detailed explanation of the symbols is given in Table 1 at the end of this description. Next, the mass conser- 
vation law or the continuity equation is given by the following equation. 

35 SO 

-jl + V • (p L g L v L +p s g s v s ) = 0 ( 9 ) 
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45 



50 
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[0043] Eq. (8) does not explicitly describe C£ itself, the solute concentration in liquid phase. Therefore, by combining 
Eqs. (5) to (9), a series of equations for equilibrium and non-equilibrium solidification type alloys are derived as follows. 

~^+v L VC„ L = V(£>„ i VC„ i ) + S (10) 

[0044] Here, the coefficients for the equilibrium solidification type alloys (denoted by suffix j) are given by Eqs. (12) to 
(15). 
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TO 



15 



A. =- 



A jMj) C j +B j.l-U) 



1- 



(\-m-A jMj) )g s ( gL+gs ) 



V-^Sl I lfl-P)g L +0 -4.0, )& X0-/?)& +g s } 



B - +b jmj) 



F-/^, +&} 



(12) 



(13) 



20 



£ _ A j,k(j)Cf +B JMJ) 
1 <\-P)g L 



(14) 



25 



£,= & 



(15) 



[0045] Likewise, the coefficients for the nonequilibrium solidification type alloys (denoted by suffix i) are given by Eqs. 
(16) to (19). 



30 



35 



40 



45 



(16) 



f 



' {l-fi)g L 



C,= ^ U ' 



^-.L.oW | B iMi) 



V 



A- MO 



i-gr-gr ] 1 -' 



(17) 



' Q-P)g L 



(18) 



50 



D,= 



1 i}-fi)g L 



(19) 



55 [0046] Furthermore, solidification contraction p is defined by 
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Ps • Pl 



(20). 



5 (3) The Relationship between Temperature and Volume Fraction of Solid 

[0047] Given the volume fraction of solid g s , the corresponding solute concentration of liquid C n can be determined 
and then the temperature is defined as a function of Cn . Thus, 



10 



T= T{d\,C L 2 ) 



(21). 



[0048] Here, it is assumed that the liquidus temperature of a multi-alloy system during solidification is determined by 
the superposition of the temperature drops in the binary phase diagrams of mother metal and each alloy element. Then, 
is the relation of Eq. (21) can be expressed by Eqs. (22) and (23) (Ref. (15)). . 



N 



1 T k-\ + lL m nMn)( C n ~ C n 9 k (n)-l) (22) 
Where , 

5-! = 3« - " + Z - qf W ) (23, 



20 



25 



The details of each symbol are shown in Table 1 . Also, N denotes the number of alloy elements. 
[0049] Next, differentiating Eq. (22) with respect to time and substituting the above-mentioned Eg. (10) to it, the tem- 
perature-volume fraction of solid relationship of Eq. (24) is obtained. 
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dT 
dt 



+ v, ■ VT= S 



(24) 



[0050] Where. S is given by Eq. (25). 



f N 



s= 



V«=i 



dt 

f N 



/V 



\tr=[ 



o-mgL s 



(25) 



[0051] A n , B n , C n and D n in the above equation are given by the aforementioned Eqs. (12) to (20). 
(4) The Darcy Equation 

55 [0052] It is well known that the flow of the interdendritic liquid is described by Darcy's equation as follows (refer to 
p.234of Ref. (14)). 



13 



BNSDOCID: <EP_ 



_0909598A1J_> 



EP 0 909 598 A1 



K 
V Sl 



(Vf + X) 



(26) 



[0053] Where, the vector v L is the flow velocity of interdendritic liquid, u is the viscosity of the liquid, K is the perme- 
ability, P is =he liquid pressure and X is the body force vector such as gravitational/centrifugal forces that includes the 
io electromagnetic body force (the Lorentz force). 

[0054] Besides, K is a constant determined by dendrite morphology (the geometrical structure of dendrite), and is 
given by the following equation of Kozney-Carman (Ref. (17)). 



75 



K = 



V-gsY 
fsl 



(27) 



20 

[0055] Here, Sb is the surface area per unit volume of the dendrite crystals (termed specific surface area), f is a 
dimensionless constant having the value of 5 as determined from flow experiments through porous media. Although K 
is basically a tensor quantity having anisotropy, it is obtained by the following two methods. 

25 Method 1 : Determination of the Permeability by A Dendrite Solidification Model 



30 



35 



[0056] In order to determine Sb in the equation of K, it is necessary to define a concrete dendrite morphology and 
take into account the solute diffusions in solid and liquid. 

[0057] Kubo and Fukusako (Ref. (18)) made a dendritic solidification model where a dendrite is modeled to comprise 
trunks and branches with cylindrical shape and tips with half-spheres as shown in Figure 5, and derived conservative 
law of a solute element at solid-liquid interface. Introducing the super-cooling phenomena (refer to pages 152 and 266 
of Ref. (14)) due to the curvature effects at the cylindrical and spherical interfaces, they derived the equation of Sb and 
showed that the calculated values of the permeability K agreed well with the measured values. 
[0058] In Figure 5, the shaded portion denotes high solute concentrated region where the solute atoms are rejected 
from the interface. Also, d is the diameter of dendrite cell, r the radius of half-sphere of dendrite tip. 
Thereupon, extending their method to the aforementioned nonlinear multi-alloy model, these inventors derived the fol- 
lowing equation. 
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Sb=a 



N 



Am* 



i 



2 ff (4iMfl{£ ^kin)) 
n=l n 



dt a LS T 



(28) 



50 



[0059] Where, a is the correction factor introduced to correct the errors of various physical properties. C n * is the sol- 
ute concentration of liquid at solid- liquid interface and can be approximated as 



55 



C L n = C L n (average sloute concentration of liquid) « C L n . 



<|> and a LS areshown in Table 1. 

[0060] From Eq. (28), Sb and then K at time t+At can be calculated from C n , gs and the solidification rate d Qs/di at 
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tme t. 

[0061 ] Furthermore, it is known from Stereology that the relationship between Sb and dendrite cell diameter d is given 
by the following equation. 

Sb = -ip (29) 

[0062] Where, <|> is the configuration factor that is $ = 1 for sphere and (J> =2/3 for cylinder (The Application of Powder 
io Theory, Maruzen Co., Ltd. (1961), p. 87, p. 132). Since the neighboring dendrite cells collide with each other when g s 
becomes about 0.7, the value of d at g s = 0.7 is calculated from Eq. (29) and is regarded as the size of dendrite cell at 
the time of the completion of solidification .i.e., dendrite arm spacing. 
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Method 2: Determination of the Permeability by An Empirical Method 

[0063] Substituting Eq. (29) into Eq. (27) and taking f=5, Eq. (30) is obtained. 



(1 " Qs) Z d Z 

f 2 (30) 

20 180<t» g s 

When the dendrite is chunky, 4» may be set 4> = 1 (Ref. (19)). The dendrite arm spacing, DAS, is determined by local 
solidification time t f as the following empirical relation (from p. 146 of Ref. (14)). 



das = A{t f ) (31) 



[0064] Where, A and n are materials constants and the diameter of dendrite cell d can be evaluated from Eq. (31) by 
30 substituting the elapsed time from the beginning of solidification instead of t f . 

[0065] Although Eq. (30) is a simplistic equation, it can not describe the accelerated solidification phenomena at the 
central portion. Also, it lacks rigidity when treating segregation. 



(5) The Momentum Equation 

[0066] The flow of the liquid in a complete liquid region is described by the Second Law of Newton, i.e., 
(mass)x(acceleration)=(force acting on body) . In other words, this is equivalent to the conservation law of momentum 
that says "the time change of the momentum (=massxvelocity) is equal to the force acting on the body" as shown in Eq. 
(32). 

^-(PV) = ZF, (32) 

d t 



The right side of Eq. (32) is the sum of pressure, viscosity force, body force, etc. Thereupon, the momentum equation 
regarding the liquid flow during the solidification process can be expressed by Eq. (33). The meanings of symbols are 
so given in Table 1. 
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^(pLg^)+V.(p i g,v I vf)=V.( A Vvf)-W+^-^^( /=] ^ 3 ) 
at K 



(33) 



Eq. (33) is solved so as to satisfy the continuity equation of Eq. (9). The suffix i denotes each component in a given coor- 
dinates system (for example, v 1 = v x , v 2 = v y v 3 = v z in (x, y, z) orthogonal coordinates system). The left side of Eq. (33) 
is the inertia term into which the volume fraction of liquid g L was introduced for the convenience when combining with 
Eq. (9) (the continuity equation). The first on the right is the viscosity force term, the second is the pressure term, the 

15 third is the sum of various body forces and the fourth is the resistant force term due to Darcy flow. 

[0067] Eq. (33) enables it possible to treat the whole region as one without distinguishing the liquid, the solid-liquid 
coexistence and the solid regions. That is, that the equation becomes a usual momentum equation by setting g L = 1 
and K= a large number, that it becomes Darcy resistance-controlled in the solid-liquid coexisting zone (inertia and vis- 
cosity forces become negligibly small) and that the liquid velocity v / becomes practically zero by setting ^ = a large 

20 number) (Ref. (20)). 

(6) The Treatment of Pearlitic Transformation 

[0068] In the case that the surface of solidifying shell is strongly cooled, the pearlitic transformation may take place 
25 because of the temperature drop at the surface layer. The volume fraction of pearlite g p is given by Eq. (34) from the 
nucleation and growth theories in the continuous cooling process. 



g p = 1 - exp(- V w ); V ex = \f{T)(t - i) z dx (34) 



[0069] Where, Vex is the extended volume of pearlite particle, t is the time and T is the temperature. The function of 
temperature f (T) is obtained from the TTT diagram for a given steel (Ref. (21 )). The latent heat of the pearlitic transfor- 
35 mation is given by pL P dg P /dt {L P : latent heat of the transformation) and is incorporated into the source term of Eq. (2). 

A-2. Discretization of The Equations 

[0070] The above equations for describing the solidification phenomena have been formulated by using the symbols 
40 of the gradients (v( ) or grad ( )), the divergences (v • ( ) or div ( )), etc. of scalars and vectors in order to make the 
operation of equations easier, to express in concise form and to be available to all coordinates systems. Next, in order 
to carry out the computation by computer, it is necessary to express these equations according to each coordinates 
system such as orthogonal and cylindrical systems and then implement volume integrals with regard to the volume ele- 
ment such as shown in Figure (7) to write them down into concrete forms. This process is called discretization. 
45 [0071] Discretization was done on the basis of the method by Patankar (Ref. (20)) in this invention. Below its outline 
is stated. 

[0072] In general, when a scalar or a vector of a physical quantity is represented by <|>, the conservation law regarding 
<t> is expressed by Eq. (35). 



^ + v.( P v^).v.(rv^) + s ( 35) 

o t 



Where, p is density, v is velocity, T is diffusion coefficient regarding <t>, S is source term regarding <|>. The velocity field 
must satisfy the condition of continuity which is given by the following Eq. (36). 
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-^ + V-(p V) = 0 (36) 

5 dt 



Eqs. (35) and (36) are expressed by defferential form. Therefore, takir-g the case of 3 D orthogonal coordinates as the 
70 example, carrying out the volume intergral iUldtdxdydz (t is time) in the volume element as swon in Figure 8 and tidying 
up with respect to <)>, a siries of Eqs. (37) to (46) are obtained (refer to p. 101 of Ref. (20)). In Figure 8, the shaded area 
denotes the control volume element and the points denoted by the circles are called grid points. Fe, Fw, Fn, Fs, Ft, Fb 
at the control volume faces e, w, n, s, t, b (t and b denote the faces parallel to the paper) denote the incoming and out- 
going of a physical quantity fy. 

15 

apt> P = Za nb 4> nb + b (37) 



Where, the suffix P denotes the defined position of the physical quantity <)> (not necessarily be geometrical center) within 
20 the volume element. The suffixes nb denote 6 neighboring definition points (E,W,N,S,T,B). These are called grid points. 
In addition, a nb {a E ,a w ,a N ,a s ,a T .a B ) are the coefficients given by Eq. (38). 

a nb =D nb A(\P nb \) + (±F nb ,o: (38) 

25 

Further, a P on the left of Eq. (37) is given by Eq. (39). 



^P=Jlci nb +af -S P AV 



(39) 



30 



old 

a P = 



old 



p =p p 



AV 
At 



(40) 



35 The source term b on the right of Eq. (37) is given by Eq. (41). 

b = S c A\/ + af X'* (41) 

40 Where, the upper suffix 'old' means the value at time t in the computational step from time t to time t + At. AV is the vol- 
ume of the volume element. D nb is the diffusion term regarding the physical quantity <|> in each face (e, w, n, s, t, b) of 
the volume element and is given by the following Eq. (42). 



^^"f^ (42) 
°nb 



T nb and A nb are the diffusion coefficients and the areas of the faces defined at these control volume faces, respectively. 
6^ correspond to the distances (6x) e ,(6x) w , • • • between the grid points as shown in Figure 8. F nb are the flow terms 
so representing the incoming and outgoing quantities of § passing through the faces (e, w, n, s, t, b) and is given by the 
following Eq. (43). 

55 

The signs of Eq. (38) are defined plus (+) when flowing into the volume element and minus (-) when flowing out. Also, 
the symbol '< )' of the second term on the right side of Eq. (38) means to take the larger of ± F nb or 0. For example, 
consider the case that <|> is taken as temperature. Fw becomes effective because of flowing in at the face w, thus Tp is 
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20 



40 



45 



influenced by the upstream side of temperature Tw. On the other hand, -F e becomes ineffective because of flowing out 
at the face e, thus T P is uninfluenced by the downstream side of temperature T E . Thus, this operation can take into 
account the physical rationality (note that the effect of the fluid flow is included in the function A(|P|) as well as men- 
tioned below). P nb is the Peclet number that describes the degree of the relative effects by the flow and diffusion, and 
is defined by the next Eq. (44). 

Pnt = (44) 
u nb 

[0073] The function A(|P|) is given by Eq. (45). 

A(\P\)= <0,(1-0.1|P|) 5 > (45) 

[0074] Considering that the source term S becomes a function of <|> in general, it is linearized as the following Eq. (46), 

S = S c + S p <t> P (46) 



in which S c and S p are the constants determined by the meaning of the equation. 

[0075] The results of discretization of the various equations described in the above section A-1 are presented at the 
end of this description. Regarding the coordinates system, the orthogonal curvilinear coordinates system was used so 
as to fit the profile of the cast piece that is elongated and curved to casting direction as shown in Figure 9. Each discre- 
25 tization equation has been written down according to this coordinates system. Since the cylindrical and Cartesian coor- 
dinates (3D orthogonal) systems are included as simple cases of the orthogonal curvilinear coordinates, the 
discretization equation can be applied to these systems as well with minimum amount of corrections, for instance by 
eliminating unnecessary terms from the equation. According to the above manipulation, each discretization equation 
becomes applicable to various cast piece profiles and cross-sections. 

30 

A-3. Analysis of The Defects 
(1) The Macrosegregation 

35 [0076] The average solute concentration of the solid-liquid coexisting zone is defined by Eq. (47) for an equilibrium 
solidification type alloy (j type), as shown in Figure 3 (b). (Note that g L + g s + g v = 1 ) 

~C j =Uc I fp L g L +C S jP s g s ) (47) 



Also it is defined by Eq. (48) for a non-equilibrium solidification type alloy (i type) from Figure 3 (c). 

C, ==(Cf p L g L +P,g,\l'crdg,) (48) 

so When ~C n > C° , segregation is defined positive; and when V n <Cn , negative. 
(2) The Influence of Dissolved Gas in Melt Steel 

[0077] It is well known that the dissolved gas in the melt steel concentrates in interdendritic liquid phase as solidrfica- 
55 tion proceeds and causes gas-caused microporosity. Here in this description, the method of analysis is described 
according to the paper of Kubo et al (Ref. (19)). 

[0078] Since the main cause of gas porosity in cast steels is CO gas, it is asssumed that CO is the only gas source. 
Then, CO gas forms by the next reaction. 
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C L + O l = CO (gas) 



(49) 



The equilibrium CO gas pressure P co is given by Eq. (50). 

5 

P co = C r OJK co (50) 

Where, C L is the carbon concentration in the liquid phase, O l is that of oxygen and Kqq is the equilibrium constant. 
io [0079] Provided that oxygen also reacts with Si that is usually added as a deoxidizer to form Si0 2 (solid) (the effect 
of Mn is neglected), the mass conservation laws regarding C and O are given by the following Eqs. (51) and (52). 

C L 9 L + C s g s +a c • P co Qv /T=c o (51) 



75 



0 L g L + O s g s + a 0 • P co g v /T + (1-y)&SiO 2 = O 0 (52) 



Where, g v is the volume fraction of gas porosity. 

The carbon and oxygen concentrations in the solid phase are given by the following equations using the equilibrium par- 
20 tition ratios. 

0 S -f<Fe-cO L (53) 
25 O s = k Fe _ 0 0 L (54) 

Similarly, Eqs. (55) to (58) hold with respect to the reaction of Si and O. 

Si L +20 L = Si0 2 (solid) (55) 

30 

K Si0 2 =S>L ' O l 2 (56) 
35 S'lQl + $ i s9s + y ASi °2 - s 'o (57) 

S' S = *Fe-S/S'L (58) 



40 [0080] P co and g v during solidification can be obtained by solving the above simultaneous equations (Eqs. (50) and 
(52) to (58)). Also, it is clear that the information about the formation of the non-metallic inclusion Si0 2 can be obtained 
as a result of computation, although it is not mentioned in this description. The meanings of symbols and the physical 
properties of the materials used in this description are listed in Table 3 at the end of this description. 

45 (3) The Effective Radius And Growth Law of Porosity 

[0081] As shown in Figure 6, the porosity is considered to form at the neck of dendrite where the local free energy 
becomes minimum (Ref. (19)). Then, defining the effective radius r of the porosity, r is modeled as follows. 
[0082] Now, it is assumed that one liquid space exists between a pair of dendrite arms and that these small spaces 
so are three-dimensionally distributed as shown in Figure 6 (b). Then, taking D as 3D mean value of the distances between 
the dendrite arms and n as the number of the liquid spaces, the volume fraction of liquid g L is approximated by Eq. (59). 



55 
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4 3 3 

3 7rrn 4*r 3 

Also from Figure 6 (c), the relationship between r, D and dendrite cell size d is shown by Eq. (60). 

2r + d = j2D (60) 

[0083] Then, from Eqs. (59) and (60). Eq. (61) of the radius r is obtained. 

]_ 

0.43865(1-^)3^ 
r-a d j- ( 6i) 

1-0.8773(1-^)3 



25 [0084] However, considering the difficulty to accurately evaluate r for the real complicated dendrite morphology, a cor- 
rection factor a d was introduced and set to be 0.7 by experience. From the above equation, it can be seen that as g s 
increases, r decreases and that as the cooling rate increases, r decreases. 

[0085] In the case that the dissolved gas is not taken into account, the equilibrium gas pressure becomes 0. Even in 
this case, the shrinkage-caused porosity form when the liquid pressure becomes less than the critical pressure. In this 
30 case, the equation regarding the growth of once formed internal porosity is given from the continuity equation of Eq. (9) 
as follows (the influence of the deformation of solid is ignored). 



dSv= — PL dg s + — V(p L g L v L ) ( 62) 
Ps Pl 



The first term on the right is the contribution due to solidification contraction, and the second term is the contribution 
40 due to the divergence of liquid phase. When dg v >0, the porosity grows; when dg v < 0, the porosity reduces (or dis- 
appears). 

A-4. Method of The Numerical Analysis 

45 [0086] All the discretization equations and various sub-equations necessary for the computation have been obtained 
as above-mentioned. There are seven equations in total that comprise the basis of the solution method: Namely, the 
energy equation, the solute redistribution equation (although there are as many numbers of equations as alloy ele- 
ments, they are counted as one for brevity), the temperature-volume fraction of solid equation, 3 components of the flow 
velocity equations regarding Darcy or momentum equation, and the pressure equation. Correspondingly, there are 

so seven major variables: Namely, the temperature T, the volume fraction of solid g s , the solute concentration of liquid C „ , 
3 conponents of the flow velocity vector and the liquid pressure. Therefore, the solution can be obtained by solving 
these discretization equations under the initial and boundary conditions. Since these variables have interaction with 
each other (it is called coupling), it is necessary to obtain the converged solution by an iterative computational method. 
[0087] Furthermore, the microscopic phenomena characterized by the permeability K determined by dendrite mor- 

55 phology, the liquid density p L as a function of solute concentrations and temperature and the formation of interdendritic 
microporosity (g v ) get deeply involved in the solidification phenomena of the heat and fluid flow in macroscopic scale. 
With regard to the solid velocity, theoretical or measured values are used. 

[0088] The numerical method developed by these inventors is described below according to the flowcharts (Figure 1 0 
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(a) and Figure 10 (b)) . 

w 

©. Set initial and boundary conditions, etc. ( Step S1 of Figure 10 (a)) 
The iterative convergency steps from time t to time t + At are as follows. 

5 ©. Calculate pressure and velocity fields of the liquid phase for a given field pattern of liquid, solid and mushy 

zones and given field patterns of permeability and liquid density (Step S2 of Figure 10 (a)). Here, either the Darcy 
equation or the momentum equation (including Darcy flow resistance) may be selected. In the case of the former 
method, solve the pressure equation to obtain the pressure field and calculate the velocity field by using thus 
obtained pressure field. In the case of the latter, the velocity and pressure fields are calculated by an extended SIM- 

io PLER method described later on. 

®. Judge the criterion of the formation of microporosity from the pressure field (step S3 of Figure 10 (a)). If pores 
form, calculate the volume fraction of the porosity and their sizes (step S4 of Figure 10 (a)), 
(g). Based on the calculated liquid flow field, the volume fraction of porosity and the heat extraction rate from the 
surface of the cast piece, solve by coupling the energy equation, the solute redistribution equation and tempera- 

15 ture-volume fraction of solid equation to obtain temperature, volume fraction of solid and solute concentration of liq- 
uid (step S5 of Figure 10 (a)). 

<S). Based on the calculated fields of the temperature, the volume fraction of solid and the solute concentration of 
liquid, calculate specific surface area Sb and the size of dendrite cell d and subsequently the permeability K using 
the dendrite solidification model (step S6 of Figure 10 (a)). 
20 (D. Calculate the liquid density based on the temperature and the solute concentration of liquid (step S7 of Figure 

10(a)). 

@. Check if the liquid pressure field is converged (step S8 of Figure 10 (a)). If converged, calculate the macroseg- 
regation from Eqs. (47) and (48) (step S9 of Figure 10 (a)): if not, go back to ©and repeat the computations. That 
is that since the permeabilities calculated in © and the liquid densities calculated in <g) affect the flow velocity field 
25 of liquid phase, the computations are repeated using these values. 

[0089] Next, the method for calculating the pressure and the flow velocity distributions using the momentum equation 
in above ©is described in detail. 

30 ©. Set the velocities at time t as the initial values (step S1 of Figure 1 0 (b)). 

©. Calculate the coefficients a p a N , a s , a T , a B , a w . Se. b of the velocity discretization equations and v 1t v 2 , v 3 
(step S2 of Figure 10 (b)). 

®. Calculate the coefficients of the pressure discretization equation (Eq. (E. 86)) (step S3 of Figure 10 (b)). 
<g). Introduce the boundary conditions for pressure (step S4 of Figure 10 (b)). 
35 (g). Calculate the liquid pressure field from the pressure discretization equation (step S5 of Figure 10 (b)). 

<g). Calculate the velocity field from the velocity discretization equation based on the calculated pressure field {step 
S6 of Figure 10(b)). 

@. Check if the velocity field satisfies the condition of continuity (step S7 of Figure 10 (b)): If not satisfied, obtain 
the corrected values of pressure by solving the pressure correction equation (Eq. (E. 1 18)) and then correct the 
40 velocity field from Eqs. (E. 1 12) to (E. 117) (step S8 of Figure 10 (b)). And return to ©. 

[0090] As above mentioned, the solution method for calculating pressure and velocity fields using the momentum 
equation was newly developed by these inventors in which various modifications/expansions were incorporated on the 
basis of SIMPLER method, one of the numerical methods in the field of heat and fluid flow analyses. Therefore, Ihis 
45 method is named the Extended SIMPLER method in the sense that the method was extended to embody the solid-liq- 
uid coexisting zone. 

[0091] Furthermore, TDM A method (Tridiagonal-matrix algorithm, p. 52 of Ref. (20)) suitable for iterative convergency 

calculation was used to solve the various discretization equations presented in this description. 

[0092] Finally, the features of the computer program of the numerical method of this invention are described below. 

50 

(1) The above-mentioned numerical method is applicable to the continuous castings with various cross-sections 
and machine profiles (vertical, vertical-bending and bending casters, etc.). Also, various analytical functions are 
available from a simple case of calculating only temperature and volume fraction of solid to the highest level of 
incorporating the effects of the deformation of cast piece, electromagnetic body force (Lorentz force), etc. in addi- 
55 tion to all aforementioned equations. Accordingly, an appropriate calculation level may be selected depending on 
the purposes; thus the highest level is not always necessary. 

The levels of the numerical analysis defined in this description are as follows. 
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Level 1 : The governing equations are the energy eq. and the Darcy eq. 



The function is porosity analysis. 

Experimentally or theoretically determined relationship between temperature and the volume fraction of 
5 solid is used. 

Level 2: The governing equations are the energy eq., the temperature-volume fraction of solid eq., the solute 
redistribution eq. and the Darcy eq. 

io The function is macrosegregation analysis. No calculation of porosity. 

The multi-alloy model is used. 



20 



Level 3: The porosity analysis is added to level 2. 



15 Level 4: The governing equations are the energy eq., the temperature-volume fraction of solid eq., the solute 

redistribution eq., and the momentum eq. 

The function is macrosegregation analysis. No porosity analysis. 

The multi -alloy model is used. The Darcy flow resistance is included in the momentum equation. 



Level 5: The porosity analysis is added to 



Furthermore, this program is installed with the functions to handle the deformation of cast piece and the elec- 
tromagnetic force. Also, the influences of the heat of Joule and the latent heat of pearlitic transformation are taken 
25 into account in the energy equation. The output includes the metallurgical informations in microscopic scale such 

as macrosegregation, microporosity. etc. in addition to temperature, volume fraction of solid, liquid pressure and 
velocity in macroscopic scale. 

(2) The above-mentioned numerical method adopts a non-steady solution method which makes it possible to ana- 
lyze through the whole process from the time of pouring into dummy bar box to the steady state and further to the 

30 final stage of solidification after the stoppage of pouring. It is also possible to analyze the effects of the changes in 
the casting speed and cooling condition, etc. throughout the process. Whether or not the steady state is reached is 
judged by observing the temperature changes. 

In conventional methods handling this kind of problem, it is common to use a steady method by the use of spa- 
tial coordinates system where the equations are written using a coordinates system fixed in space and the steady 

35 state solution is obtained by iterative computation (thus, the computational domain is fixed in space). However, 
there are the shortcomings in them that the important part of non-steady state can not be analyzed. On the con- 
trary, the above-mentioned non-steady method possesses the advantages that enables it possible to accurately 
respond to the changes in various external factors (thermal, mechanical, etc.). 

(3) With respect to a vertical-bending type, etc., the cast piece undergoes bending deformation. Accordingly, the 
40 topologies (the distance, area, volume etc.) of the object for analysis changes as shown in Figure 9 (b), and also, 

the components of gravitational force changes as shown in Figure 9 (a). To correspond to these changes, they are 
updated as time proceeds. 

(4) The boundary condition at the surface of cast piece is given either by the heat transfer coefficient h at the sur- 
face (thereafter called h-method) or by the surface temperature Tb itself (thereafter called Tb-method). With h- 

45 method, Tb response is obtained; with Tb-method, h response is obtained. For example, when a particular surface 
temperature distribution is desirable, obtain h by using Tb-method, and determine the corresponding cooling con- 
dition from the relationship between h and the cooling condition (such as the quantity of water-spray). 

(5) In the mushy region where the liquid pressure drop occurs, the flow direction of the liquid phase can be regarded 
mainly as one-dimensional flow towards the casting direction. Therefore, solving Darcy Eq. (26) with respect to the 

so body force X z in Z direction (casting dir.), Eq. (63) results. 

x ' m -^*nr v ^ (63) 

55 

Accordingly, after the pressure and velocity field is obtained (with no porosity formation), define the P distribu- 
tion optionally to prevent the formation of porosity (for example, by taking the pressure gradient from the position of 
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P=0 to the crate end as dPIdZ = 0), and calculate X 2 (the sum of the gravitational force in Z dir. and the Lorentz 
force) from Eq. (63). Thus, it is possible to obtain the required electromagnetic body force distribution (Lorentz 
force). 

(6) Considerable amount of input data are given by external functions. For example, the operating conditions (cast- 
5 ing temperature, casting speed, surface cooling condition, etc.) are given by the functions of time, position, etc. 

(7) The merit of the nonlinear multi-alloy model developed in this invention is to be able to expand the applicability 
of the above-mentioned numerical method by fitting to the nonlinearity in the phase diagrams. Thus, the method 
can be applied to many important commercial alloys regardless of ferrous, nonferrous, stainless steel, etc. For 
example, as for the carbon steels of C=0.1 to 0.51% with peritectic reaction, the relationship between temperature 

io and the volume fraction of solid can be obtained by neglecting the peritectic reaction and smoothly approximating 
the 8 and y solidus lines. It is of course possible to apply to the carbon steels less than 0.1 %C. 

[0093] The items of the above (1) to (6) are not included in Figure 10 because it becomes very complicated. Further- 
more, in this computer program, the solid phase in the solid-liquid coexisting zone is assumed not to flow (yet the defor- 
75 mations accompanying the bending/unbending and the soft-reduction are acceptable). Regarding this assumption, 
even if the solid crystals are assumed to move in the region with quite a low fraction of solid (to say about 0.3 or less), 
the resulting effects can be neglected. This is shown by the best modes for carrying out the invention (mentioned later) 
where the liquid pressure drop due to the interdendritic liquid flow is very small indeed. Thus, the above assumption is 
adequate. 

20 

A-5. Computational Example of The Numerical Analysis 

[0094] As an computational example, the steel of 0.72% C-0.57% Si-0.70% Mn-0.02% P-0.01% S-remainder Fe 
(wt%) was chosen which has a marked tendency that the liquid density decreases as the solute concentrations of the 
25 interdendritic liquid increases during solidification. The numerical analysis was done for the solidification process of the 
steel cast into the mold of 1 m diameter x 3m height as shown in Figure 1 1 (a). The initial temperature was set at 1475°C 
(superheat 13°C). The physical properties used are those of the 0.55wt% carbon steel given in Table 2 and Table 3 at 
the end of this description. 

[0095] The computation was started from the time that the mold was filled with the melt. During about 1 0 minutes until 
30 the substantial solidification from the mold wall begins, the liquid basically flows descent at the mold wall side and 
ascent at the central region and is turbulent flow. In other words, ® The flow velocity is about 10 cm/s at the rapid flow 
region, © The temperature inversion layer appears where the temperature at the central part becomes lower than that 
at theside. Thus, the flow pattern is that of turbulent, the temperature distribution becomes quickly uniform (the temper- 
ature difference is less than 2°C) and most of the superheat is lost. After that, such situation continues until about 2 
35 hours-when the liquid phase disappears and all area becomes a solid-liquid coexisting zone. In the meantime, the flow 
velocity gradually becomes small. 

[0096] Solidification begins from the bottom, spreads to the side and finally ends at somewhat upper position from the 
central part of the ingot (solidification time is 20.9 hrs). 

[0097] It is seen that the isolines of temperature and volume fraction of solid after 1 1 .5 hours are largely curved as 
40 shown in Figure 1 1 (b) and (c). This can not be expressed by mere temperature calculation, thus reflecting the effect of 
liquid flow in the solid-liquid coexisting zone which follows next. In these diagrams, the symbol C denotes shrinkage cav- 
ity, M solid-liquid coexisting zone (mushy zone) and S solid zone. 

[0098] As shown in Figure 1 1 (d) and (e), the interdendritic liquid flow is such that the flow pattern becomes descend 
at the central portion and ascend at the side portion. Although the temperature at central portion is higher than that at 

45 the outside (accordingly light), and the interdendritic solute concentrations of liquid are lower at the central portion than 
at the outside portion (accordingly heavy), the balanced liquid density at the central portion becomes heavier than that 
at outside; therefore resulting in the above-mentioned flow pattern. This flow pattern continues to the latter half of the 
solidification. As pointed out from the solidification theory by Flemings et al (p.244 to 252 of Ref. (14)), the positive seg- 
regation takes place by the flow from lower temperature portion (higher solute concentrations) to higher temperature 

so portion (lower solute concentrations): The negative segregation takes place by the reversed flow (from higher to lower 
temperature). Shown in Figure 11 (f) and (g) are the segregation distributions for C and P, respectively. The other ele- 
ments (Si, Mn, S) shows the same trend; thus reflect well the macrosegregation found in large steel ingots. 
[0099] Since the number of elements used is relatively small (7 in radius dir. x 30 in height dir. evenly partitioned), it 
is not possible to express the locally concentrated V segregation themselves at the center or A segregation themselves 

55 at the upper part of the ingot (refer for example p.244 of Ref. (14)). However, the formation processes of the segregation 
are well grasped, thus showing the validity of this simulation. 

[0100] The above calculation is the results of the strictest analysis that uses a prescribed dendrite solidification model 
and applies the momentum equation to the whole ingot without distinguishing the liquid and mushy zones (analysis 
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level 4 but no porosity analysis done). The computational error was evaluated by the difference between the heat 
extraction from the ingot Qout and the heat loss of ingot Qlost by, |(Qout-Qlost)/Qoutx 100|(%) . In the case of only tem- 
perature calculation, the total amount of errors till the end of solidification was less than 0.1%. 

5 B. The Mechanism of the Formation of the Internal Defects 

[0101] Many literatures have been published about the internal defects of cast steels. 

Figure 43 schematically shows the central defects (internal defects) that takes place in an elongated cast steel bar. The 
regions A and C in the diagram are sound due to the leeding effects of interdendritic liquid. The region B is unsound 
io because of the difficulty of the liquid feeding and results in the formation of the interdendritic microporosity along the 
centerline. It is known that this microporosity usually exhibits the form of V characters pointing the feeding direction as 
shown in Figure 43 and, in many cases, accompany V form of macrosegregation (so-called V segregation) (for example 
Ref. (34)). 

[0102] There are few papers that clearly distinguished V form of microporosity and segregation in the past literatures. 
15 For example, Pellini (Ref. (35)) has called them centerline shrinkage without distinguishing. The central defects (internal 
defects) in continuous castings o1 steels are essentially the same as those of the above-mentioned cast steels. There- 
fore, summing up the V forms of porosity and segregation as one regardless of the existence or the degree of these 
defects, they are called the central defects in this description. 

[0103] The central defects form in the case that the interdendritic liquid feeding is insufficient. Thus, it can be said that 
20 the flow of the liquid phase in solid-liquid coexisting zone (or mushy zone) plays a decisive role to the formation of the 
central defects. As the driving force that causes this liquid flow, the following factors are pointed out: 

(1) The flow due to the solidification contraction that is induced by the density difference of the solid and liquid dur- 
ing solidification. In addition, the effects of the thermal contraction associated with the temperature drops of the 

25 solid and liquid phases are included. 

(2) The flow due to the density difference within the liquid phase (natural convection). The liquid density p L depends 
not only on the temperature but also on the solute concentrations in the liquid phase, as shown in the following Eq. 
(64). 

30 Pl-Pl<Ci.c£.- • -.C^.T) (64) 

(3) Forced flow due to mechanical deformation from outside such as bulging, unbending, soft-reduction, etc. This 
is more comprehensible when imaging the flow of the water upon squeezing or bending the water absorbed 

35 sponge. Additionally, intensively cooling cast piece to cause the thermal contraction enters to this classification. 

[0104] These inventors conducted a series of preliminary numerical analyses to investigate the above-mentioned fac- 
tors (1) and (2) on the central defects. The results are summarized as follows. 

40 (J). The macrosegregation forms prominently in large steel ingot. This is because the interdendritic liquid flow 
occurs in the wide range for a long period. When decreasing the size of the ingot for the same alloy, the range of 
the mushy zone becomes narrow, but the flow pattern shows the same tendency as those in Figure 1 1 (d) and (e). 
However, since the solidification time is shortened, this flow is limited to relatively small range and the segregation 
does not form practically. This coincides with the experiences. Thus, when the solidification rate is increased, the 

45 natural convection caused segregation due to the density difference in liquid phase barely takes place. 

@. In continuous casting, the flow pattern is not that of the natural convection type due to the difference in liquid 
density as seen in the examples described later on, but that of the simple solidification contraction caused flow in 
the casting direction. With slabs, the cooling intensity often differs in lateral direction. Even in this case, as far as 
"norma! solidification" takes place, the degree of the difference in macrosegregation within the plate is relatively 

so small; thus, acceptable from practical point of view. This also is attributed to the rapid solidification rate. The Darcy 
flow pattern in the normal solidification is, as shown in Figure 1 2 (a), slightly outward (in the figure, the outward flow 
is somewhat magnified to emphasize this pattern). Besides, in the vicinity of the central portion where the central 
defects form, the flow velocity in the casting direction is overwhelmingly larger in comparison with that in the thick- 
ness direction; thus, the flow velocity in the thickness direction is negligibly small. 

55 

[0105] The above-mentioned numerical analyses were done for the whole cast piece ranging from the meniscus to 
the final solidification position using Level 3 analysis. From a viewpoint of the whole Darcy flow pattern, the influence of 
the outlet flow from the nozzle is small. 
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[0106] From above, it can be said that the major internal defects taking place in continuous castings are the V pattern 
defects that form in the final solidification position within the cross-section and that the solidification contraction caused 
flow is most deeply involved as a governing factor. 

[0107] Next, the mechanism of the formation of the V pattern defects is explained. 
5 [0108] Since the main flow in the mushy zone, elongated along the casting direction, occurs in the casting direction, 
most of the liquid pressure drop due to the Darcy flow takes place in that direction. In particular, the pressure drop 
becomes largest at the central portion of the cross-section and its neighborhood. Then, if the liquid pressure P reaches 
to the critical condition given by Eq. (65), porosity forms (p.239 of Ref. (14)). 



Where, Pgas is the equilibrium partial gas pressure within the porosity in equilibrium with the dissolved gas in the liquid, 
is a LG is the surface tension at the liquid-porosity interface, and r is the radius of curvature of the porosity given by Eq. 
(61). The porosity is arranged along the V patterns as shown in Figure 12 (b). In the case that the V segregation devel- 
ops, it is considered that triggered by the formation of the porosity the Darcy flow in the casting direction shifts from the 
normal pattern of Figure 12 (a) to the flow pattern of as indicated in Figure 12 (b). In other words, the liquid flows in 
along the V patterned voids from the lower temperature portion (higher solute concentrations) to the higher temperature 
20 portion (lower solute concentrations) and results in the local positive segregation band where the solute concentrations 
are higher than the average; thus forming the V segregation. It is considered that once the porosity is formed, such liq- 
uid flow takes place in a simultaneous manner with the formation of porosity. 

[0109] If the flow velocity from the lower to higher temperature is increased and reaches to the condition given by Eq. 
(66), the local remelting phenomenon would occur (p.249 of Ref. (14)). 



25 



30 



Vt • Vt 

1+ — <0 (66) 

9t / at 



[0110] When such remelting is brought about, the Darcy flow resistance of the remelted portion becomes lower than 
that in the surroundings and the flow is progressively increased resulting in more remelting. Consequently, the V seg- 
35 regation would appear more notably. The degree of the segregation depends on the extent of this channeling phenom- 
enon (according to the second term on the left of Eq. (66)). 

[01 1 1 ] In order to examine the above argument of the formation of the central defects, a carbon steel was melt by a 
high frequency induction furnace and cast into the tapered dry sand mold of 32 to 30mm diam. x 350mm long as shown 
in Figure 44. Further, as a means to enhance the interdendritic liquid feeding, the mold was placed in a pressurized 

40 cylindrical vessel as shown in Figure 44, and pressurized by argon gas after pouring. 

[0112] The chemical composition of the cast specimen is shown in Table 4. The casting temperature was 1560 to 
1580°C, the pouring time was about 10 seconds. The oxygen and nitrogen contents were 50 to 120 ppm. As to the No.1 
specimen cast in the atmosphere, the thermocouples were inserted at 3 locations along the center as shown Figure 44 
to measure the temperature changes during solidification. The measured temperatures were shown in Figure 45. 

45 [0113] In addition to the casting experiments, the numerical analyses by this invention were carried out from the 
beginning of pouring to the end of solidification and compared with the experiments by tracking the formation process 
of internal defects. The physical properties used are shown in Tables 2 and 3. The chemical compositions were set 
those of Table 4. As to the mold, the thermal conductivity was set 0.0036 cal/cms°C, the specific heat 0.257 cal/g°C, the 
density 1 .5 g/cm3. As to the insulating material of the riser, the thermal conductivity was set 0.0003 cal/cms°C, the spe- 

50 crfic heat 0.26 cal/g°C the density 0.35 g/cm3. 
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Table 4 





The chemical compositions of the pressurized 


5 




cast specimens ( wt%) 






No 


C 


Si 


Mn 


P 


S 




1 


0.40 


0.28 


0.42 


0.024 


0.018 


10 


2,3 


0.31 


0.18 


0.22 


0.045 


0.017 



The calculated temperatures (denoted by the broken lines in Figure 45) at each location of No.1 specimen have well 
agreed with the measured values. 

The macrostructure of No.1 specimen etched by 0.4% nital solution is shown in Figure 46. Figure 46 (a) shows the 
15 schematically sketched V patterns by naked eye observation and Figure 46 (b) a part of as etched macrostructure 
where the dark etched central V defects are clearly visible. The macrostructure consists of the columnar structure at 
very near surface and the fine equiaxed structure. Figure 46 (c) shows the position of microstructure and the measuring 
position of Vickers hardness. The microstructure and the result of the Vickers hardness measurement are shown in Fig- 
ure 47 and Figure 48, respectively. In Figure 47, the needlelike white portion is ferrite and the dark etched matrix is 
20 pearlite. 

The dark etched band flowing from upper left to down right of Figure 47 shows that the ferrite is scarce at the band and 
so the carbon concentration there is higher than that in the surroundings. Measuring the Vickers hardness across this 
flow as shown in Figure 46 (c), the resulting hardness was higher in the V band with increased pearlite than that in the 
surroundings as shown in Figure 48. In addition, the hardness once decreased at the vicinity of the V band and subse- 
ts quently increased to rise the right as shown in Figure 48. This is considered to be attributed that when the V band is 
formed, the higher solute concentration liquid flows in (in this case, from left side) along the V band. 
[01 1 4] Also, it was confirmed by the color check inspection on the central cross-section that the microporosity distrib- 
uted along the V pattern. The volume of the shrinkage cavity (of Figure 46 (a)) was about 1% of the total volume of the 
casting. This is quite smaller than the solidification contraction 4% of the casting, showing that most of the defects exist 
30 as the microporosity in the V band. 

[0115] It was thus confirmed that the V patterned defects consist of the microporosity arranged in V characters and 
V segregation bands (the positive segregation). 

[0116] Shown in Figure 49 is the results of Level 3 analysis for the No.1 specimen where the formation process of 
microporosity was analyzed. The porosity distribution after solidified (Figure 49 (b)) is in good agreement with the real 
35 V pattern (Figure 46 (a)). From the numerical results, the internal porosity form 55 seconds from the start of pouring 
and the distribution of volume fraction of solid at that time is shown in Figure 49 (a). 

The range of the porosity is denoted by the hatched area. By Level 2 analysis assuming no porosity formation, the pres- 
sure drop due to the Darcy flow became the biggest 63 seconds after from the start of pouring at the position of 75mm 
from the bottom: the pressure was -20.5 atm. Considering the above results, the computations were performed by 
40 changing the atmospheric pressure in the range from 10 to 25 atm. The results indicated that the critical pressure for 
the formation of porosity is 20 atm. Thus, as the pressure was increased, the volume fraction of porosity decreased. 
Therefore, it is expected that the porosity completely disappears with more than 20 atm. 

[0117] Based on the above investigations, in No.2 specimen, the pressurization was started (the volume fraction of 
solid at the center is about 0.3) and kept at 10 atm from after 30 seconds from pouring to the end of solidification. The 
45 macrostructure is shown in Figure 50. Although the volume fraction of porosity was decreased compared to No.1 spec- 
imen, the V defects are prominently visible. 

[0118] On the other hand, as shown in Figure 51 of the macrostructure of No.3 specimen cast at 22 atm, the sound 
region without V segregation and porosity extended from 30mm to 130mm showing that the pressurization worked 
effectively. 

so [0119] With these specimens No. 2 and No. 3, Level 3 numerical analyses were conducted with the results shown in 
Figure 52 (the chemical compositions are by Table 4). The porosity decreases to some extent at 10 atm and disappears 
at 20 atm. The internal defects formed below the riser (Figure 51). This is because the amount of the melt was less and 
the shrinkage was deepened. Also, in the numerical analysis, the formation of shrinkage cavity at the riser was not 
treated strictly (to strictly treat this, the partitioning of the computational elements around the riser should be fine 

55 enough, but for this particular case, not done because of the problem concerning the display of the results). 

[0120] From above, it is clear that the internal defects can be eliminated by pressurization, reconfirming experimental 
results published in the past. However, in these past experiments, the effect of pressurization was not studied theoret- 
ically and quantitatively; therefore the experiments were inadequate. For example, in Ref. (34), the authors expressed 
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a negative perspective about the effect of pressurization in practical use, quoting that the central segregation appeared 
more prominently in the pressurized casting at riser. In this reference, the sizes of the casting (3 inch rectangular cross- 
section x 24 inch long), chemical compositions, casting temperatures, pressurization conditions at riser, measured data 
of temperature during solidification and the results of the observations of internal defects are all given; accordingly, it is 

s possible to compare with the numerical analysis. 

[0121] Thereupon, these inventors conducted Level 3 three-dimensional analyses of this cast steel with the results 
shown in Figure 53. It was found that the effect of pressurization is small at 4.2 atm, and that at least 20 atm is required 
to eliminate the central defects. In this connection, in terms of the criterion Eq. (69) of the formation of porosity, the rela- 
tionship between the liquid pressure drop and the formation of porosity is schematically illustrated in Figure 54. In the 

io figure, porosity forms at the critical volume fraction of solid gs*. From the above verification, it is rational to think that the 
reason why the central segregation contrarily appeared more prominently in the above Ref. (34) is because when the 
pressurization effect on the riser is insufficient, the porosity formed and then triggered by the formation of the porosity, 
the high solute concentration liquid in the vicinity of the porosity flowed in. In any event, aside from the detailed discus- 
sion on dendritic scale, it can be said that if the feeding effect is sufficient, the central defects do not form. 

15 [0122] In conclusion, the internal defects form in the region where the liquid pressure in solid-liquid coexisting zone 
(mushy zone) becomes less than the critical pressure during solidification, and it is possible by the use of the numerical 
method by this invention on the basis of the solidification theory to calculate the position of the defects in continuous 
castings. 

[0123] Furthermore, it can be said that "When the microprosities are suppressed, the macrosegregation is sup- 
20 pressed simultaneously." Thus, it is important to eliminate the microporosity, or to say more precisely, not to give a 
chance of the formation. To accomplish this, it is necessary to restrain to the minimum the liquid pressure drop associ- 
ated with the Darcy flow in the casting direction in the vicinity of the central region of thickness (the final solidification 
zone) and to hold it more than the critical pressure defined by Eq. (65). 

25 C. Calculation of The Electromagnetic Force 

[0124] Various methods are conceivable to exert the electromagnetic body force (Lorentz force). For example: Method 
of applying DC magnetic field and DC current, method of utilizing the distant propulsion force by linear motor type, etc. 
Thus, an appropriate method is available depending on the shape of the cross-section of cast piece, the exerting posi- 
30 tion, the magnitude of required force, cost of the equipment, etc. 

[0125] . Here, the calculation method is described about the former case. As shown in Figure 2, the electromagnetic 
body force (Lorentz force) f acting in casting direction is given as an outer product of the current density vector J in lat- 
eral direction and the magnetic flux density vector B in thickness direction by Eq. (67). 

35 

f = J X B (67) 



40 J is expressed from Ohm law as Eq. (68). 



J= a E = -crV^ (68) 

45 

Where, E is the electric field strength, v4» is electric potential gradient, a is electric conductivity. And, the electric poten- 
tial distribution is obtained by the following Eq. (69) (from Eq. (3.4) of p.31 of Ref. (22) and the above Eq. (68)). 

50 V • (aV<|>) = 0 (69) 

4>, is obtained by solving Eq. (69) using the boundary condition of the electrical potential defined at electrodes. Iron is 
nonmagnetic above the Curie point (about 770°C) and can be regarded approximately same as that of the air. There- 
55 fore, it is relatively easy to exert an uniform static magnetic field to the solid-liquid coexisting zone. By incorporating the 
calculated f into the body force term X in the Darcy equation of Eq. (26) or the momentum equation of Eq. (33), and 
performing numerical analysis, its efficiency can be evaluated. 

[0126] Further, the heat of Joule Qj is given by Eq. (70) and is taken into account in Eq. (2). 
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Qj = cjJ 2 (J/m 3 s) = 0.238889 x 10" 6 aJ 2 (cal/cm 3 s) 

[0127] Best modes for carrying out the invention by this invention is described below. 

5 

A. On the Vertical-Bending Continuous Casting of Round Billet 

[0128] As shown in Figure 13, the continuous casting machine of the best mode No.1 for carrying out the invention 
comprises water-cooled copper mold 5, tundish 3, nozzle 4 and the electromagnetic booster 1 to exert electromagnetic 
io body force onto the solid-liquid coexisting zone (mushy zone) of the cast piece. 

As shown in Figure 23, the electromagnetic booster 1 is the apparatus to generate the electromagnetic body force 
toward casting direction which comprises superconductive coils or electromagnet to generate DC magnetic field and 
electrodes to flow DC current. 

[0129] 1%C-1%Cr bearing steel with 300mm diameter was chosen. Because the bearing unit receives repeated load 
is under high speed, excellent fatigue and wear resistances are required. Accordingly, among many specialty steels, bear- 
ing steel is the one that the most severe qualities are required about the cleanliness, the uniformity of structure, etc. This 
steel has a wide solidification temperature range, is easy to bring about the central segregation resulting in the forma- 
tion of coarsened carbides and causes the quality deterioration of the low service life. The chemical compositions were 
set 1%C, 1%Cr, 0.2%Si, 0.5%Mn, 0.1%Ni, 0.01%P and 0.01%S. The physical properties used for computation are 
20 given in Table 2 and Table 3 at the end of this description, the linearized data of Fe-C phase diagram is shown in Figure 
14. The relationship between the temperature and the volume fraction of solid using the nonlinear multi-alloy model of 
this steel is shown in Figure 15 (a). 

[0130] When the volume fraction of solid g s approaches 1 in the nonlinear multi-alloy model, the coefficients A n , etc. 
(Eqs. (12) to (19)) of the solute redistribution Eq. (10) becomes infinity To avoid this inconvenience upon computation, 

25 the solidification was assumed to be completed when g s =0.95 (Ref. (16)). In this connection, the latent heat of fusion 
was corrected so as to release 100% at g s =0.95. Hence, the latent heat of fusion was assumed to evenly evolve, con- 
sidering that the mushy zone is elongated towards the casting direction. That is, the value of 68.4 (cal/g) obtained by 
dividing the heat of fusion 65 by the fraction solid 0.95 was regarded as the apparent latent heat of fusion. The dissolved 
oxygen condenses in interdendritic liquid as the solidification proceeds, giving rise to the change of the equilibrium CO 

30 gas pressure as shown in Figure 16 (a) (refer to Eqs. (49) to (58). Assumed that there is no CO gas bubble). The phys- 
ical properties used are given in Table 3. From the above figure, it is obvious that Pco falls off with the decrease of O 
content. 

[0131] The computational domain were divided 10 evenly in radius direction (Ar = 1 .75cm), and the partition length of 
the elements in casting direction was set Az= 5cm. Prior to the computations, the number of elements in radial direction 
35 was examined where the temperature variation is steep. It was found that no essential difference was observed with 
more than 8 elements. The same examination was done with respect to the casting direction: Thus the number of par- 
titions was determined as above mentioned. 

[0132] With respect to the correction factor a of the specific surface area of dendrite Sb (Eq. (28)), a was determined 
a = 1.2 so as to coincide with the measured values of the dendrite arm spacing (DAS) of 1C-1.5Cr steel of Eq. (71). 

40 

DAS = 523 (average cooling rate, °C/min) -0.55 (*im) (71) 

[0133] The above equation was obtained using the average cooling rate during solidification temperature range of 
45 (1453-1327=126°C). The temperature of the melt steel flowing in from the tundish was set constant and the radiation 
heat from the meniscus surface was ignored. 

[0134] The liquid flow pattern in the upper part of the melt pool becomes complicated because of the outlet flow from 
nozzle, convection flow due to the temperature difference within the pool. Thus, the flow is basically turbulent and the 
temperature difference within the melt pool becomes small. Also, as already mentioned, the behavior of the liquid pres- 

50 sure drop in the mushy zone elongated in casting direction is most crucial. From this point of view, the effect of the liquid 
flow within the bulk liquid pool is exceedingly small. Accordingly, if we focus on the problems of internal defects, the flow 
within the melt pool does not necessarily need to be analyzed in detail. Considering these points, the solution method 
by the Darcy equation was used instead of solving the momentum equations that requires an excessive computation 
time. According to the Darcy solution, the fluid flow inside the bulk liquid pool becomes modest; as a result, the thermal 

55 diffusion due to the convection becomes small. To compensate this, the thermal conductivities of the liquid pool and the 
mushy region with the volume fraction of solid less than 0.05 was apparently multiplied by 5 that of the liquid. (This 
method is frequently used to calculate the temperature of continuous castings (for example, refer to Ref. (24)). However, 
in these computations, the flow analysis is not done and the errors brought about by neglecting the Darcy flow is cor- 
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rected by using apparently increased thermal conductivity.) As an example, the comparison is shown in Figure 17 
between the case that only temperature was calculated and the case that the above Darcy solution method was used. 
It can be seen that compared to (1) of only temperature calculation, the solidification at the center begins earlier and 
the mushy zone is elongated by the influence of the flow-in of higher temperature liquid from upstream. From this exam- 
ple, it is obvious that the Darcy flow analysis is necessary even in a macroscopic scale. 

Hence, the computational results regarding the conventional operating conditions are shown in No.1 and No.2 of Table 
5, Figures 18 and 19. The computations were done applying Level 2 and Level 3 analyses. 



70 



Table 5 



Best mode No.1 for carrying out the invention: The analytical results of 1C-1Cr bearing steel of vertical continuous 

casting 


(Casting speed 0.6m/min, casting temperature 1473°C (superheat=20°C)) 


No./ Casting 
method 


Computational 
condition: Poros- 
ity analysis 


M length (m) 


Z (max) (m) 


Pmax (atm) 


Porosity g v (%) 


1 Conventional 


Not done 


16.05 


20.9 


-8.6 




2 Conventional 


Done 


15.45 


20.3 


-0.10 


9.5%ati=1 
5.8% at i=2 


3 Eprocess: 
Lorentz force of 
8G exerted in 
range 19.5 to 
21m from menis- 
cus 


Done 


16.05 


21.0 


-0.03 


No porosity 


Note 1 : M length is from the position of the volume fraction of solid g s =0.01 to the crater end at g s =0.95. 
Note 2: Zmax is the length from the meniscus to the crater end. The position that liquid phase disappears (in other 
words g s + gv =1) is defined the crater end and the next element where the liquid phase exist is defined the crater 
end element. 

Note 3: Pmax is the liquid pressure at the crater end element. 

Note 4: Symbol i in the volume fraction of porosity g v denotes the element number in radial direction. 
Note 5: The segregation forms in the case that the porosity forms. 
Note 6: E process is the method by this invention. 
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30 



35 



40 



45 



[0135] The heat transfer coefficient h at the mold surface was gradually changed from 0.02 to 0.01 (cal/cm2s°C) to 
prevent the breakout (Figure 18 (b)). 

[0136] In the secondary cooling zone by water-mist spraying, the surface temperature of the solidifying shell was set 
uniformly to 1 125°C. Then the h can be obtained as a response. The heat flux from the surface is given by the product 
of h and the difference between the surface temperature and the ambient temperature. The boundary condition was 
changed from the mist cooling to natural radiation cooling at the position where the cooling intensity by radiation 
becomes larger than that by the mist cooling (refer to Figure 18 (b) and (d)). 

[0137] The liquid pressure at the crater end element (the farthest from the meniscus, the final solidification position) 
becomes -8.6 atm in the Level 2 analysis of the computation No. 1 . However, such a negative pressure can not be real- 
ized in practice, thus porosity forms according to the following critical condition for the formation of porosity. 



50 



P < P 



2a 



LG 



CO 



(65)(Rewritten) 



[0138] Pco increases to the maximum value of 0.9 atm with increased volume fraction of solid. On the other hand, the 
55 term -2a LG /r increases approximately to the negative value of -1 .2 atm for this particular case. Accordingly, unless the 
liquid pressure (in the side of the higher volume fraction of solid with bigger pressure drop) becomes less than the pres- 
sure of P (absolute pressure) = 0.9-1 .2 = -0.3 atm, the porosity does not form. 

[0139] In the computation No.2 by Level 3 analysis with the porosity formation taken into account, the liquid pressure 
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became less than this critical value resulting in the formation of porosity. The relationship among the liquid pressure P, 
the gas pressure Pco and the volume fraction of porosity g v after the formation of porosity is automatically adjusted to 
satisfy Eq. (65) and the already mentioned Eqs. (49) to (58). In this computer program, the Darcy flow is allowed even 
in the situation that the porosity exists. In this way, 5 to 10% porosity formed in the central range of about 6cm (20% of 
5 the diameter). 

[0140] g v in here represents the mean value of the volume fraction of porosity in a considerably larger volume element 
in comparison with the order of the dendrite arm spacing. Similarly with respect to the central segregation, the com- 
puted values are also those of the mean values in the volume elements. Hence, the segregation is not brought about 
upon computation. However, this does not mean that the V segregation doe not form, but the V segregation does locally 
10 form in fact as already mentioned. 
Some comments are given below. 

1) As shown in Figure 18 (a), P increases approximately linearly in low fractions of solid region (say less than 0.2). 
Thus, the pressure drop is very small, which means that even if the starting point of the mushy zone at upstream is 

is changed to some extent, its effect on the pressure drop at the higher volume fractions of solid region associated 
with the porosity formation is almost negligible. From this, it is understood that the strict flow analysis by the 
momentum equation in the melt steel pool is not necessary. The distributions of the solid, the liquid and the mushy 
are shown in Figure 20. 

2) The Darcy flow is descent with the maximum value of -2.8mm/s and becomes slower as it goes to upstream 
20 because the width of the stream becomes wider (similar to the flow of river). In the upper part of the liquid pool, 

ascent flow is observed which is the natural convection that results from the higher temperature at the central part 
compared to the side. The variation of the body force X (gravitational force) is shown in Figure 18 (c). The body 
force X becomes smaller at the crater end side. This is attributed that the effect of the condensation of the lighter 
solute elements (all except for Ni) than the liquid Fe is greater than the effect of temperature drop, thus resulting in 

25 smaller liquid density p L . As already mentioned, the driving force to cause the Darcy flow is the contraction associ- 
ated with solidification, and the flow is downward almost uniformly as shown in Figure 12 (a). The flow pattern in 
the vicinity of the crater end for the computation No.1 is shown in Figure 21. The flow channel becomes narrower 
toward the crater end and the flow velocity in the radial direction gradually becomes smaller compared to that in the 
casting direction (in the vicinity of the crater end, the flow in the radial direction is practically negligible). 

30 3) The degree of the segregation of the alloy elements at the central portion is within the computational error of a 
few percents, i.e., practically no segregation (however, note that the V segregation takes place as above men- 
tioned). 

4) The permeability K regarding the Darcy flow is one of the important factors when evaluating liquid pressure drop. 
In this description, the two methods for determining K were described. The dendrite arm spacing obtained by these 

35 methods are shown in Figure 22. In the curve (a) obtained by Eqs. (28) and (29), DAS is smallest at the surface, 
becomes larger as it goes inward but on the contrary becomes smaller at the central portion. This is attributed to 
the accelerated solidification at the final stage of solidification, as is seen that the shell thickens rapidly as it 
approaches the crater end (Figure 18 (b)). On the other hand, as shown in the curve (b) by Eq. (31), the local solid- 
ification time t f becomes the biggest at the center and hence, the dendrite arm spacing DAS becomes the biggest 

40 at the center as well. The accelerated solidification phenomenon appears at the last stage of solidification more 
clearly in prescribed large steel ingots. It is also found in continuous casting as reported in Ref. (25), thus it is a 
common phenomenon. With respect to the distribution of DAS in thickness direction, it is reported that it becomes 
small conversely at the central part of the continuous casting of 6063 aluminum alloys (the diameter 203mm, cast- 
ing speed 0.1m/min) of Ref. (26). 

45 

[0141] From the above, it is obvious that the solution method by this invention theoretically evaluating d and K by the 
use of Eqs. (27) to (29) reflects the solidification phenomena more strictly. This is one of the reasons to use these equa- 
tions. [The above reversible phenomenon can not be grasped by Eqs. (30) and (31). This is because the history of solid- 
ification rate dQs/di is not considered in these equations.] 
so [0142] Next, the results of the computation No. 3 by Level 3 analysis for the case that applied the electromagnetic 
force by this invention are shown in No. 3 of Table 5 and in Figure 24. The conceptual schematic diagram applying this 
method to the vertical continuous casting of round billet is as shown in Figure 23. The Lorentz force is given by Eq. (72) 
as the product of uniform DC magnetic flux density B x in X direction and the DC current density J y in y direction that 
flows in the central part of the solid-liquid coexisting zone. 

55 

f 2 = - J y B x (72) 
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[0143] In the light of the P distribution in Figure 18 (a) and the required Lorentz force obtained from Eq. (63), the 
Lorentz force of f z = -54900 (dyn / cm 3 ) (8 times of gravitational force, 8G) was exerted onto the range from the upstream 
vicinity of the position of P becoming 0 to the crater end, i.e., the range of 19.5 to 21 .0m from the meniscus. As a result, 
the liquid pressure drop nearby the crater end with high volume fraction of solid is relaxed to hold a positive pressure of 
about 1 absolute atm and the porosity did not take place as can be seen from Figure 24. In other words, the continuous 
casting without the internal defects can be produced by exerting the Lorentz force larger than the above value. 
[0144] With respect to the soft -reduction method by the prescribed references, it has been interpreted such that the 
soft-reduction method tries to suppress the interdendritic liquid flow toward casting direction by giving onto the solidify- 
ing shell the reduction gradient corresponding to the solidification contraction and thereby tries to reduce the central 
defects. However, this can be reinterpreted such that it relaxes the liquid pressure drop occurring in the casting direc- 
tion. In this sense, since it is possible to reduce the required Lorentz force by concurrently applying the soft-reduction, 
it is effective to attach soft -reduction gradient by placing rolls between the round billet and rigid frame la as shown in 
Figure 2 (d). This will be discussed in detail later on. 

B. On The Vertical- Bending Continuous Casting of Thick Slab 

[0145] The vertical-bending continuous casting of thick slab is described as the best mode No.2 for carrying out the 
invention. 

[0146] Since the central defects in thick slab, for example, high grade thick steel plate used in ocean structures orig- 
inate cracks and thus cause the quality deterioration, it has been studied energetically as an important problem that 
influences the quality. The central segregation appears more prominently in higher carbon content steels with a wide 
solidification temperature range. Thereupon, 0.55% carbon steel of JIS S55C (AISI 1055) was chosen. The chemical 
composition was set 0.55%C, 0.2%Si, 0.75%Mn, 0.02%P, 0.01%S. The relationship between the temperature and the 
volume fraction of solid obtained by the nonlinear multi-alloy model is shown in Figure 15 (b) and the physical properties 
in Tables 2 and 3 and in Figure 25. Si is used as a deoxidizer. The oxygen content was set 0.003 wt%. The outline of 
the caster is shown in Figure 26 and the operating conditions given in Table 6. 



Table 6 



The specification and the operating conditions of the vertical -bending caster 
used for the best mode No.2 for carrying out the invention 


Mold length 


1.2 m 


Length of vertical section (including mold) 


3 m 


Bending radius of curvature 


8 m 


Dimensions of slab 


220 mm thick x 1500 mm width 


Casting speed 


1 m/min 


Superheat of melt steel 


15°C 


Oxygen content in melt steel 


0.003 wt % 



[0147] The slab undergoes plastic deformation when it passes through bending and unbending zones. Assuming a 
simple bending mode and regarding the position of neutral axis as unchanged considering that the radius of curvature 
is large enough compared to the slab thickness, the strain in casting direction z z becomes maximum at the surface with 
the value of e^l 10/8000=1.375%. The radii of curvature shown in Figure 26 were set so that the total bending strain 
of 1 .375% was obtained by gradually bending with 5 steps by about 0.275% per each step. They were set similarly with 
respect to the unbending zone. In this computer program, assuming a simple bending /unbending deformation mode 
as above mentioned and regarding the cast piece as a complete plasticity body (elastic strain is ignored), the effect of 
the plastic deformation is included as solid deformation velocity in various governing equations. 
[0148] The computational domain was partitioned uniformly into 1 9 elements throughout in thickness direction, con- 
sidering the non-symmetric nature by bending (partition length Ax = 22cm / 19). The partition length in casting direction 
was set Az = 10cm. Since the width of slab in lateral direction is considerably larger compared to the thickness, two- 
dimensional analyses were performed. The correction factor a of the specific surface area of dendrite Sb (Eq. (28)) was 
set 1 (no correction). 

[0149] First, the results of Levels 2 and 3 analyses for the conventional operating conditions are presented in Table 7 
and in Figures 27 to 31 . The computations No.1 and No.2 are those for the conventional methods at present. The heat 
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transfer coefficients at surface h (cal / cm 2 s°C) are set as: 

h = 0.03-0.001 5V2 in the mold 
h = 0.015 for Z = 1 to 3 m 
5 h = 0.010 for Z ^ 3 m 

Where, Z is the distance from the meniscus. The surface temperature and solidified shell thickness changes are as 
shown in Figure 27 (d) and (b), respectively. 



10 

Table 7 





The best mode No.2 for carrying out the invention: Analytical results for 0.55% carbon steel of vertical -bending 

caster 


75 


(Casting speed is Im/min, casting temperature is 1500°C (superheat=16°C)) 




No./ Casting process 


Computational con- 
dition: Porosity 

analw 
cti idiy. 


M length (m) 


Z (max) (m) 


Pmax (atm) 


Porosity g v (%) 


20 


1 Conventional 


Not done 


8.7 


18.6 


-4.7 






2 Conventional 


Done 


8.5 


18.4 


-0.3 


8% at center ele- 
ment; Diam of pore : 


25 


3 Eprocess 












Range:18.0 to 18.6m 

Magnetic flux density: 
0.7(T) 


Done 


9.1 


19.0 


-0.1 


No porosity 


30 


Current density: 
1.47x10 6 (A/m 2 ) 

Lorentz force: 15G 
(G: gravity) 












35 


3 Eprocess 














Range:same as 
above 


Done 


9.4 


19.3 


0.78 


No porosity 


40 


Magnetic flux density: 
0.5(J) 

Current density: 
2.058 x 10 6 (A/m 2 ) 












45 


Lorentz force: same 
as above 












Note: The meanings of Symbols, etc. are the same as in the note of Table 5 



[0150] When the volume fraction of solid g s becomes more than 0.6 in the computation No.1 by Level 2 analysis, the 
so liquid pressure drops sharply and becomes a negative pressure of -4.7atm at the crater end. This is attributed that the 
permeability K decreases rapidly as shown in Figure 27 (c). The casting directional component X of the gravity 
becomes 0 at the range more than Z= 16m, hence there is no feeding effect by the gravitational force (Figure 27 (c)). 
Consequently, the porosity forms. 

[0151] In the computation No.2 with porosity analysis (refer to Table 7 and Figure 28), about 8 vd.% of porosity were 
55 formed in the range of 1 1 mm (5.2% of the thickness) about the center. The size of the porosity is about 50 um in diam. 
That means that a severe V segregation takes place that accompanies the porosity at the central region. Attention must 
be paid that the liquid pressure distribution differs in the computation No.1 and in the computation No.2 with porosity 
formation taken into account. The large negative pressure of the computation No.1 is not able to take place in reality, 
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and the real pressure distribution becomes as shown in Figure 28 as a result of the porosity formation. 
[0152] Next, described are the computations No.3 and No.4 where the electromagnetic force by this invention was 
applied. In reference to the information about the required Lorentz force distribution (by Eq. (63)) obtained along with 
the computation No.1 , the parameters were set as below in the range of more than Z= 18m from the meniscus where 
5 the liquid pressure drops significantly: 

Range Z = 18.0m to 18.6m from meniscus 

DC magnetic flux density in thickness direction of slab B = 0.7 (T) 

DC current density in width direction of slab 

10 

J = 1.47 x 10 6 (A/m 2 ) 
To generate, 

15 Lorentz force toward casting direction of slab 
f = J x B = 1.029 x 10 6 (N/m 3 ) 
(equivalent to 1 5G, 1 5 times the gravity) 

20 

[0153] For this, the potential difference at both ends of the computational domain in the width direction was set as 
follows. 

E = J x 0.01 /a = 1 .47 x 10 6 x 0.01 / 7.0 x 10 5 = 0.021 (V) 

25 

[01 54] The electric conductivity a is taken the average in mushy zone (Figure 32). In this example, the electromagnetic 
booster is installed in the horizontal zone of Figure 26. It is also possible to extend the applied range and thereby reduce 
the required Lorentz force. 

[0155] Thus obtained results are presented in Table 7, No. 3 and in Figures 29 and 30. It is obvious from Figure 29 

30 (a) that the liquid pressure drop was relaxed to -0.1 1 atm (absolute pressure of 0.89 atm) at the crater end element, 
thus no porosity forms. The central segregation is a few percent within the level of computational error and thus essen- 
tially does not exist. The whole solidification profile and the Darcy flow pattern in the vicinity of the crater end are shown 
in Figure 30. The flow pattern is normal in the range where the electromagnetic force was exerted and in the unbending 
zone. In the unbending zone, the deformation mode is that of tensile at the free side (inside the curvature) and that of 

35 compressive at the fixed side (outside the curvature) about the center axis. Therefore, as a result, the interdendritic liq- 
uid is squeezed out by the reduction in thickness at the free side (this is reversed at the fixed side) and the liquid flows 
from the free to the fixed side (this phenomenon was clarified by the preliminary computations, but is omitted for want 
of space). However, in this example, such liquid flow due to unbending was not observed. This is because the bending 
strain is basically small with the maximum of e 2max = 1 4% at tne surface, and further becomes smaller at the central 

40 region. From the above, it can be said that the unbending deformation does not influence the macrosegregation when 
the deformation is approximated by the simple bending deformation mode. Heat of Joule was admitted to generate to 
some extent in Lorentz force applied zone as can be seen from Figure 29 (c). This leads to the elongation of Zmax 
(length from meniscus to crater end, often called metallurgical length) from 18.6m to 19.0m (40cm elongated). 
[0156] When the product of current density J and magnetic flux density B is constant, the resulting electromagnetic 

45 body force (Lorentz force) f becomes constant too. However, it is desirable to make J as small as possible and increase 
B upon operation, because the vicinity of the center of slab remelts by the heat of Joule when J is too large. Here, on 
the contrary, B was decreased down to 0.5 (Tesla) and J increased up to 2.058 x 10 6 (A/m 2 ) to generate the same 
Lorentz force, and the effect of heat of Joule was investigated. The results are shown in Table 7, No.4 and in Figure 31 . 
Compared to the computation No.3, the effect of the heat of Joule becomes further larger, and Zmax is elongated 70cm 

so from 18.6m to 19.3m. At the crater end element, the liquid pressure is held at the positive value of 0.78 atm. Thus, it is 
understood that this level of the heat generation has no problem. However, if the remelting of central portion occurs, it 
takes time to re-solidify and the mushy zone is elongated again and the pressure drop occurs again. In such a case, it 
becomes meaningless to apply Lorentz force. 

[0157] From this, it is desirable to increase the magnetic flux density and lower the current density. The apparatus 
55 using super conductive magnet that is able to produce a high magnetic flux density would be more advantageous from 
a viewpoint of economy, the save of space, etc. compared to conventional magnet. This will be mentioned again later. 
Also, in this example, the DC current was supplied through the whole thickness at the sides of the slab. But, practically, 
it is sufficient to supply only in the vicinity of the central region where the Lorentz force is required and thereby make it 
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possible to reduce the heat evolution by the heat of Joule. 

[0158] It is clear from the above example that the internal defects can be eliminated by exerting the electromagnetic 
force by this invention for thick slabs as well. 

5 C. On the Vertical -Bending High Speed Continuous Casting for Thick Slab 

[0159] High-speed casting is taken up as the best mode No.3 for carrying out the invention. In general, the productivity 
(given by output tons per caster per month) is determined by non-operating time, preparation time, dimensions in cross- 
section, casting speed, etc. Among these, the important factors are the dimensions in cross-section and the casting 

10 speed that are closely associated with the quality: Enlarging the cross-section is not worthy from a metallurgical point 
of view, thus much efforts have been paid to raise the casting speed. Thereupon, the case is described below that 
applies this invention to a high speed casting of slab that has increasingly been preferred. The specification and the 
operating conditions are such that the casting speed was set 2 m/min and all other parameters except for the cooling 
condition were set the same as those in the best mode No.2 (of Table 6) for carrying out the invention. 

15 [0160] The computational results by conventional method are presented in Table 8, No.1 and in Figure 33. 



Table 8 



20 


The best mode No.3 for carrying out the invention: Analysis results of 0.55% carbon steel of vertical-bending high 

speed casting 




(Casting speed 2 m/min, casting temperature 1500°C (superheat=16°C)) 


25 


I'Uw V— /cl o Lilly |ji uucoo 


vUl 1 IfJUlCtMUl Idl OUIKJI 

tion: Porosity analy. 


ivi lengin \it\) 




Pmax (atm) 


Porosity g v (%) 




1 Conventional 


Not done 


14.5 


33.1 


39.2 




30 


2 Conventional 


Done 


12.6 


31.2 


1.15 


1 5% at center, pore 
diam 65 ^m; 

5% at both sides of 
the center, pore 
diam 60 ^m 


35 


3 E process: 














Lorentz force of 1 5G 
for 2=30.2 to 31.7m: 
B=1.33(T) 
J=7.775x10 5 (A/m 2 ) 


Done 


14.8 


33.4 


0.16 


No porosity 


40 


Lorentz force of 34G 
for Z=31.7to 33m: 
B=3.0(T) 

J=7.775x10 5 (A/m 2 ) 












45 


4 E process: 














Lorentz force of 8G 
for Z=30.8to 33.1m: 
B=1.38(T) 
J=4x10 5 (A/m 2 ) 


Done 


14.6 


33.2 


5.1 


No porosity 


50 


Reduction grad.: 
0.1 mm/30.8 to 
33.1m 












55 


Note: The meanings of Symbols, etc. are the same as in the note of Table 5 



[0161] Zmax, 33.1m (crater end length or termed metallurgical length), becomes 1 .8 times longer in comparison with 
Table 7, No.1 and the liquid pressure drop was increased to the negative value of -39.2 atm. In the computation No.2 



34 

BNSDOCID: <EP 0909598A1 J_> 



EP 0 909 598 A1 



(Table 8, No.2 and Figure 34) with the porosity analysis done, about 5% (35mm range about the center, 16% of thick- 
ness) to 1 5% (at the center element) porosity were formed. The size of the porosity is similarly in the range of about 60 
urn to 65 jim. The Lorentz force equivalent to 22G on the average is required in the range of Z=30.2m to 33.1 m to elim- 
inate the porosity. Hence, dividing the negative pressure region into two zones, the Lorentz force was exerted as follows 
5 (Level 3 analysis). 

No.1 zone: Lorentz force equivalent to 15G is exerted in the range 2= 30.2 to 31 .7 m from meniscus. For this, the 
parameters are set as follows. 

w DC magnetic flux density B = 1 .33 (T) 

DC current density J = 7.775 x 10 5 (A/m 2 ) 
Potential difference for the analytical domain in slab's lateral dir. 
E = Jx0.01/cr = 0.0111 (V) 

75 No.2 zone: Lorentz force equivalent to 34G is exerted in the range Z= 31 .7 to 33.1 m from meniscus. For this, the 

parameters are set as follows. 

B = 3.0 (T) (increased) 
J = 7.775 x 10 5 (A/m 2 ) (unchanged) 
20 E = 0.01 1 1 ( V ) (unchanged) 

[0162] The results are shown in Table 8-No. 3 and in Figure 35. The liquid pressure is held at P=-0.16 (atm) (absolute 
positive pressure of 0.84 atm) at the crater end: Therefore, no porosity or V segregation occurs. The metallurgical 
length Zmax was elongated from 33.1m to 33.4m. This is because the solidification was a little delayed by the influence 
25 of the heat of Joule. In this case, the average Lorentz force of 22G was applied over 2.8m toward casting direction. How- 
ever, it is desirable to reduce the applied range as well as the Lorentz force from the viewpoint of economy or equip- 
ment. 

[0163] Thereupon, as a next logical step, it was tried to relax the liquid pressure drop by using the soft-reduction as 
a supplementary means and thereby reduce the required Lorentz force (by utilizing mechanical or magnetic attractive 
30 force, see Figure 2 (d)). 

[0164] As a preparation, computations were conducted to examine the effect of soft-reduction with the results shown 
in Figure 36 (a) to (c). Figure 36 (a) shows the reduction distribution necessary to completely compensate the solidifi- 
cation contraction. Upon the calculation ol the reduction distribution 8 , the position where the volume fraction of solid 
g s at the center element becomes 0.1 (or any number would do) was set as the reference position (in this particular 
35 case, Z=25m), and the volume contraction due to solidification was calculated in the mushy zone. Finally, the 8 distribu- 
tion was obtained by taking 8 =0 at the reference position and by equating this volume contraction to 8 in the thickness 
direction of the slab. Thus, taking AS as the increment of the reduction quantity during At, 



AS = f2>(/)0 S (''M/) (73). 

i 



Where, S is the cross-sectional area perpendicular to thickness direction of volume element, p is the solidification con- 
45 traction, g s is the solidification rate, V is the volume of element and the suffix i denotes the mushy element in the thick- 
ness direction. Next, in reference to the calculated reduction distribution, the actual reduction gradients were given as 
shown in Figure 36 (b) and the computational results obtained was presented in Figure 36 (c) which shows the degree 
of relaxation of liquid pressure drop in the vicinity of the crater end. In the light of the above results, the reduction gra- 
dient of 0.10mm / 30.8 to 33.1m was given simultaneously with the Lorentz force of 8G in the same range. The result 
so is shown in Table 8-No. 4 and in Figure 37. The defects disappeared completely. The applied range of Lorentz force in 
the casting direction was shortened 50cm compared to the computation No.3 with only Lorentz force exerted, and the 
required Lorentz force decreased down to about one third, thus indicating that even slight reduction gradient is consid- 
erably effective. 

[0165] In relation to this best mode for carrying out the invention, the advantages and cares, etc. are discussed below 
55 upon applying this invention. 

(1) On the soft-reduction gradient. 

The reduction gradient given to the computation No. 4 is smaller than the value to compensate the net solidif i- 
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cation contraction, and the contraction of the cast piece due to the temperature drop toward casting direction and 
the deformation due to thermal stress are not considered. Accordingly, the real reduction quantity will be bigger 
than the value of this example. As mentioned in the BACKGROUND OF THE INVENTION, because the reduction 
gradient employed by present soft-reduction method aims to completely compensate the solidification contraction, 
s it is generally larger than that described in this description. Therefore, a possibility is pointed out in Ref. (27) that 

when the strains in the mushy zone becomes larger than a certain limit, the dendrite crystals are destroyed 
mechanically and the high solute concentration liquid is sucked to give rise to internal cracking (there are much 
questions about the detailed mechanism). 

The soft-reduction by the definition of this description is such that "the reduction quantity is fairly small (there- 
to fore less than the above-mentioned strain limit), and utmost is used as a supplementary means to relax the liquid 
pressure drop". In other words, the interdendritic liquid feeding by the Lorentz force plays a major role. Accordingly, 
it may be said that there is no possibility of the internal cracking that often brings about a problem in conventional 
soft-reduction method. 

(2) On the relation between the hydrogen induced cracking and the central segregation under sour gas environ- 
15 ment. 

Since large diameter of line pipes to transport petroleum and natural gas are served under severe environments 
such as in the underground, the sea bed, cdd district, etc., excellent properties are demanded for toughness and 
fracture characteristics as well as strength. If the hydrogen coming out from the humid H 2 S atmosphere in these 
crude oil and gas penetrates into the pipes and is trapped by the central defects (that formed during continuous 

20 casting and remained in the final production), so called HIC (hydrogen induced cracking) occurs. The corrosion 

resistance for this H 2 S is in general called sour resistance and has become an important technical matter since the 
accident by HIC of the transportation pipe in the Arabian Gulf sea bed in 1972 (Ref. (28)). 

One of the measures taken for the HIC at present is to adjust the chemical compositions to eliminate HIC, 
admitting the central segregation (and porosity) as an unavoidable reality. For example, in Ref. (29), the chemical 

25 compositions are adjusted so that the HIC sensitivity parameter P HIC given by Eq. (74) becomes less than 6 with 

special attention paid to the HIC sensitive elements C. Mn and R 

p hic - O eq + 2P* < 0.6 (wf%) (74) 

30 

Where, C* q is the equivalent carbon content given by Eq. (75). P* is segregated P concentration. S M denotes the 
degree of segregation of element M (>1). 

35 Coq=$c * C + S M n ' M n /6 + (S Cu • Cu + S Ni . A//)/15 + (S Cr • Cr+S Mo • Mo + S v • V)IS (75) 

Take the case of API (American Petroleum Institute) standard X65 class (meaning the proof stress of more than 
65000psi (448Mpa)) for example. To satisfy the specification of this steel, C and P were set significantly small values of 
C=0.03 and P=0.004 (wt%) according to this criterion, respectively and besides the compositions of Cu, Ni, etc. were 
40 severely controlled. Furthermore, special care was paid to thejhermo-mechanical treatment. 

The sensitivity parameter P HIC becomes 0.53 and Ceq =0.33 from the same reference. If there is no segrega- 
tion, these values becomes P Hlc =0.298, C eq =0.29. 

Applying this criterion to the above best mode No.2 and No.3 for carrying out the invention of 0.55% C steel 
and assuming no segregation by the use of this invention, P H | C becomes 0.71 5. P Hlc decreases to 0.365 when only 
45 C is decreased to 0.2% (however, infinitesimal additions of Cu, Ni, Cr, Mo, V are not included). This means that if 

the segregation is eliminated, the severe control over the chemical compositions, etc., become unnecessary and 
that the degree of freedom for balancing the compositions expands substantially. 

Considering that the demand for the strength of these line pipes is growing from X70 class to X80 class and further 
beyond these and at the same time stronger resistance in HIC and SSC (sulfide stress cracking), weldability and 

so so on are increasingly demanded, the meaning of expanding the degree of freedom for the balance of compositions 
is significant. The relationship between compositions and mechanical properties is omitted here for want of space. 
However, it is relatively easy to fulfill these demands, considering the fact that high strength materials have been 
developed one after another. In conclusion, it can be said that such severe demands can be fulfilled by eliminating 
the central defects via this invention. At the same time, adjustment over the compositions becomes possible; as for 

55 this case lowering C content. 

(3) With respect to the function f(T) for the 0.55%C steel in Eq. (34), it was determined from the TTT diagram of 
Ref. (30) as follows. 



36 

BNSDOCID: <EP 090959BA1 J_> 



EP 0 909 598 A1 



f(T) = 3.547 xlO" 12 • 



fr-30oY 4 * 53 f760-rY 362 



V 



100 J 



V 100 J 



(76) 



10 



15 



20 



25 



The TTT diagram obtained from Eqs. (34) and (76) is shown in Figure 38 along with the experimental data. Both 
agree relatively well. In the case of this computation, the surface temperature fell off to 540°C and resulted in 100% 
pearlitic transformation at the surface elements (thickness 1 1 .6mm) in the range: 2 = 18.7m to 22.5m (the crater 
end). The recalescence of the surface temperature Ts of Figure 33 (d) is attributed to the latent heat of pearlitic 
transformation. 

(4) Here, the attractive force generating between two coils is examined when using superconductive magnet. The 
model used is shown in Figure 39 (a). For the convenience of the computation, the coils are assumed circle and 
the total current I (= current in a superconductive wire x the number of turns) in the coil are assumed a point current 
as shown in the figure (in practice, it has a finite cross-sectional area). The cast piece exists between the two coils, 
but is regarded same as air. The magnetic flux density B 2 at the center axis of 2=b/2 is given by the following equa- 
tion (refer to standard text book, for example, Saburo Adachi "Electromagnetic Theory" , Shokohdo (first edition, 
1989), p. 79 and p.89). 



(a 2 +b 2 /4)i 



(Tesla) 



(77) 



30 



Where, [1o=4tcx10' 7 (H/m) is magnetic permeability of vacuum. On the other hand, the force in 2 direction that the cur- 
rent of coil 2 receives by the magnetic field that coil 1 makes is given by the following equation. 



F z = - ZnalB r (N) 



(78) 



35 



Where, B r is the component in r direction of the magnetic flux density on coil 2 and is given by the following equation by 
using vector potential Ae (the 9 direction component). 



40 



dz 



(79) 



45 



50 



55 



[0166] (Regarding Ae, refer to Naohei Yamada et al. "Exercise on Electromagnetic Theory" (1970), p.159 [Corona 
Company Ltd.] for example.) The results obtained using Eqs. (77) to (79) is shown in Figure 39 (b). a was fixed at 0.8m 
and B z was set 1, 2 and 3 (T). The figure shows the relationship between the pressure P (the value of F z divided by 
area ua 2 ) between the coils and the distance between the two coils. 

[0167] The aforementioned calculation uses the parameters that hypothesized actual operation, and shows that it is 
possible to control the pressure exerted on the cast piece by controlling the magnetic flux density, i.e., the coil current 
and the distance between the coils. The strength of dendritic skeleton in the mushy zone is in the range of several 
Kg/cm 2 to 50 Kg/cm 2 (p. 72 of Ref. (27)). Hence, it can be said that it is possible to give a very small reduction gradient 
by utilizing the gravity (attractive force) between the coils (refer to Figure 2 (d)). For example in the case of best mode 
No.3 for carrying out the invention, the volume fraction of solid g s at the central portion is 0.65 or more in the soft-reduc- 
tion range of 30.8m to 33.1m. Then, judging from the strength of dendrite skeleton as above-mentioned, it is possible 
to give a prescribed soft-reduction gradient by setting the distance between the coils of 0.6m and B=1 to 2 (T). Prior to 
a practical application, it is necessary to obtain in advance the empirical relationship between the magnetic attractive 
force and the reduction gradient on a real machine with the electromagnetic booster equipped (refer to Figure 40). 
Then, the magnetic attractive force may be applied for prescribed reduction gradient in reference to thus obtained rela- 
tionship. Since non-defect is guaranteed by the liquid feeding by Lorentz force with the soft-reduction used as a supple- 
mentary means to relax the liquid pressure drop, it is enough to control the magnetic attractive force within a certain 
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degree of allowance (strict control is not necessary). 



D. On the Bending Type Continuous Casting of Bloom 

5 [0168] The bending type continuous casting of bloom is taken up as the last best mode for carrying out the invention. 
The material used is the same 0.55% carbon steel as the best mode No.2 for carrying out the invention, and also the 
chemical compositions and the oxygen content were set the same values. The cross-section is that of rectangular with 
thickness 300mm x width 400mm, the bending radius of the curvature of the machine was set 15m, the length of the 
mold 1 .2m and the length of the water-mist spraying zone below the mold 4 m. 

10 [01 69] With respect to the bending type caster, the radius of curvature at the mold is the same 1 5m. Accordingly, the 
cast piece undergoes only unbending deformation at the unbending zone and the radii of curvature between unbending 
rolls were set as shown in Figure 55 such that the cast piece undergoes evenly with 4 steps of bending strain (total 
strain of 1 50mm/l 4850mm (radius of curvature at neutral axis) =0.0101). The casting temperature was set 1500°C, 
same as the best mode No.2 for carrying out the invention. The casting speed was set 1m/min. Above specification and 

15 the operating conditions are those generally used in this kind of bloom castings. 

[0170] 3-dimensional analysis by Level 3 was conducted considering that the heat flow pattern becomes that of 3 
dimensional in bloom. The computational domain was partitioned uniformly into 15 elements in radial direction (partition 
width = total thickness of 300mm/15=20mm) and the partition length in casting direction was set 150mm. Considering 
symmetric nature in width direction, the computational domain was taken half the width and partitioned uniformly into 5 

20 elements (200mm/5=40mm). The heat transfer coefficient at mold-cast piece boundary in the mold was set as 

h = 0.03 - 0.00146 V? (Z distance from the meniscus) (cal/cm 2 s°C), 
h = 0.015 in water cooling zone, 
h = 0.005 in natural cooling zone. 

25 

[0171] The physical properties used are the same as best mode No.2 for carrying out the invention. The correction 
factor regarding the specific surface area of dendrite Sb was also set <x=1 (no correction). The results of the numerical 
analysis by Level 3 for a conventional casting method are shown in Figure 56. The length of mushy zone is 14.1m, the 
crater end length Zmax is 27.9m and 5.6 vol% of porosity with the pore size of about 54 urn was formed at Z=27.82m 
30 in the center element (thickness 20mm x width 40mm): Thus it was judged that the central defects were formed. 

[0172] The Level 2 analysis was done to obtain the Lorentz force to eliminate the central defects, from which it was 
found that the Lorentz force equivalent to 18G was necessary in the range of Z=27.6 to 28.05m. Hence, the Lorentz 
force was exerted in the range of Z=27.3 to 28.05m as follows. 

35 F = J x B = 10 6 (A/m 2 )x 1.2 (T) = 1.2 x 10 6 (N/m 3 ) 

The cross-sectional area at the electrode side of the bloom was set at width 1 40mm x length 750mm, considering that 
the solid-liquid coexisting zone exists in relatively narrow range of the cross-section. The current lines become consid- 

40 erably uniform in the relatively narrow range of the central mushy zone about the center between the electrodes 
attached at both sides. The current pattern spreads in the thickness and longitudinal directions of the cast piece to 
some extent. The current that flows through the same cross-sectional area at the central part of the cast piece as that 
at the electrode was 65% of the total current (in 3 dimensional computation of the current field, all surfaces except for 
the electrodes are assumed insulated]. The results obtained are shown in Figure 57. The liquid pressure is held at suf- 

45 ficiently large positive pressure in the vicinity of the crater end: Thus, the central defects do not form. 

E. Specific Example of Electromagnetic Booster 

[0173] In this section, more detailed mechanisms of the electromagnetic booster are discussed that exerts the elec- 
50 tromagnetic force generated by DC current and DC magnetic field as already described in the above four best modes 
for carrying out the invention. Also, the specific mechanisms of the combined system of the electromagnetic force and 
the soft-reduction method will be shown. And, the mechanism to reduce the tensile force produced in cast piece by the 
electromagnetic force is described. 

[0174] First, superconductive coils are used as the means for generating the DC magnetic field, and a single pair or 
55 plural number of pairs of coils are arranged in such a way that the cast piece is placed between the coils. As for blooms 
and billets whose lengths in short side and long side of the cross-section are relatively close, the racetrack type coils 
elongated in casting direction are basically used. On the other hand, racetrack type coils broadened in lateral direction 
are used for wide slabs correspondingly. Because the superconductive coils need to be cooled to liquid helium temper- 
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ature (4.2K) at present, they are enclosed in the cooling container consisting of liquid helium, etc. Also, the reaction 
force corresponding to the Lorentz force in casting direction is exerted on the coils. Therefore, these points should be 
taken into consideration upon designing. Also, because the attractive force is generated between the coils when DC 
current flows, it is necessary to enclose the coils into highly rigid frames which are fixed by plural number of supporting 
5 columns. 

[0175] Next, as a means for supplying DC current at both sides of the cast piece, the plural number of sliding elec- 
trodes fixed at space are arranged to contact with the side surfaces of the cast piece. Thin oxidized layer consisting of 
Fe as main component forms at the casting surface. Since this layer is an insulator, it is desirable to remove it by means 
of cutting, etc. Also, as a means for enhancing the contact at the electrode-side surface boundary, the plane cutting is 
io used in this invention. Further, to prevent the re-oxidization of the cut plane, the cut surface is blocked from the atmos- 
phere by using inert gas such as argon, N 2 or reductive gas, etc. 

[0176] The soft-reduction gradient is given through a plural number of rolls, and the pressurizing devices by means 
of the fluid such as oil are adopted to the bearing unit of each roll so that an optional reduction distribution can be given 
by controlling them independently. It is necessary to reduce the distance between superconductive coils as small as 

75 possible to obtain a strong magnetic field. Thus, it is important to make the diameter of rolls smaller. As to the rolls for 
blooms and billets, it is desirable to use the rolls whose central part swelled out to be able to effectively give the reduc- 
tion onto the central mushy zone and to avoid cracking caused by unnecessary plastic deformation at the corners. As 
to wide slabs, conventional flat roll can be used. Also, it is good to use so-called divided roll which is divided into sub- 
rolls in the longitudinal direction of the roll to minimize the bending due to reduction force or thermal stress. 

20 If the Lorentz force toward the casting direction is too strong, it gives rise to the tensile stress in the cast piece with 
mushy zone and may result in internal cracking. As a means to reduce the excessive tensile force, the drawing resist- 
ance created as a result of soft- reduction may be taken advantage of. Besides that, driving device may be equipped 
with to these rolls. 

[0177] The above is the main mechanical means by which it becomes possible to control the current density distribu- 
25 tion and electromagnetic force distribution appropriately. It also becomes possible to give desired reduction gradient. 
[0178] It was shown in the above best mode No.3 for carrying out the invention that the required Lorentz force to elim- 
inate the central defects can be reduced by giving the soft-reduction gradient as a supplemental means. This principle 
can of course be applied to blooms, etc. In other words, by suitably balancing both of them, it is possible to adjust the 
balance of the drawing resistance due to the soft-reduction and the Lorentz force toward casting direction as well as to 
30 eliminate the central defects. The balance between them changes depending on machine profile, and operating param- 
eters such as casting speed, shape of cross-section of cast piece, steel grade and so on. 

In the case that both forces are balanced, the tensile force that occurs in the solid part of the solidifying shell due to the 
electromagnetic force is counterbalanced (from macroscopic point of view). When the drawing resistance is large 
enough compared to the electromagnetic force, the driving force by the rolls may be given toward the casting direction. 
35 On the contrary, when the electromagnetic force is too large, giving rise to a large tensile force in the solidifying shell, 
the rotating speed of the rolls are regulated to correspond to the casting speed. In this case, the reversed torque is 
exerted onto the rolls and works as a brake: Thus, the tensile force in the solidifying shell can be counterbalanced. In 
summary, the electromagnetic apparatus by this 
invention has the following three functions. 

40 

Function I Electromagnetic force alone 

Function n Combination of electromagnetic force + soft-reduction 

Function ru Combination of electromagnetic force + soft-reduction + roll-drive 

45 [0179] By utilizing these functions properly, individual purposes (sound cast piece with no defects / high speed cast- 
ing) can be realized. 

Specific Example 1 : Application to bloom and billet The specific example applied to steel bloom or billet is shown 
in Figure 58. With respect to machine profile, vertical-bending type or bending type is most widely used as sche- 

50 matically shown in Figure 1 . Figure 58 shows the case that the electromagnetic booster is located in the upstream 
vicinity of the final solidification zone (crater end) at the horizontal zone of cast piece. Figure 58 (a) is the cross- 
sectional plan and (b) is AA cross-sectional plan in longitudinal direction. The arrow in the diagram denotes the 
casting direction. The view from the top, BB cross-section, is shown in Figure 59. 

Symbol 6 in the figure denotes the cast piece and Symbol 102 denotes the electrode located at both sides of 

55 the cast piece that contacts with the side surface. The electrode is fixed to the frame 107 (details not shown) with 
spring 1 06 and slides against the moving cast piece. The plural number of electrodes are arranged over the Lorentz 
force zone as shown in Figure 58 (b). Each electrode is independent. The smaller the interval between each elec- 
trode, the better. 
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The connection type of the electrodes is shown in Figure 60. Figure 60 (a) is a parallel type and the current 
density that flows through each electrode is roughly equal (contact resistance is assumed constant). Figure 60 (b) 
is a series type that is suitable for the case that the current density is changed, for example, when increased density 
is favored in the downstream side. 
5 Figure 60 (c) is a mixed type of (a) and (b), and the current is set for each parallel group. Also, it is possible to 

change individual current density at electrode by changing the material of electrode. These can be selected 
depending on the situation. 

The individual electrode is stored in box 105 of insulation nature and is connected to L-type bus bar 104 and 
plate bus bar 103. The bus bar 103 shown in the BB cross-section of Figure 59 corresponds to the parallel type of 

10 Figure 60 (a). Figure 62 shows the situation where gas shield box 109 to prevent the oxidization at the electrodes 
and plane cutting milling machine 108 are attached. Figure 62 (a) is that of side view, (b) is that from the top. Sym- 
bol 110 denotes the electrode box room, Symbol 111 the milling machine box room and both rooms are divided. 
112 and 1 1 3 are gas inlets. The oxidation-preventing gas is released little by little from the gap 1 1 6 after the air in 
both rooms are exchanged with it. Plural number of cutting tools 1 14 are attached to the milling machine. 1 15 is 

15 discharging outlet of cut chips. Gas inlet box is not shown to avoid the complexity in Figure 58. 

Symbol 120 in Figures 58 and 59 is racetrack type coil wound by superconductive wires and is built in to the 
stiff frame 122. 121 is a cooling chamber for the coil and cooled to liquid helium temperature (4.2K). Since upper 
frame 117 and lower frame 1 1 8 are heated by the radiation from high temperature cast piece, etc., it is desirable to 
insert cooling chamber 123 between these frames and rigid frame 122. 

20 The upper frame 117 and lower frame 118 are supported by pillars 119, movable up and down, and can be 

locked at a specified position. Since these frames and pillars are burdened with the magnetic attractive force 
between the coils, the reaction force by reduction roll, etc., it is necessary to take the modulus of section large 
enough to reduce elastic deformation due to bending, etc. to minimum extent. Also, nonmagnetic materials such as 
stainless steel may be used. Since the rigid frame 128 is burdened with the reaction force corresponding to the 

25 Lorentz force acting in casting direction, the drawing resistant force due to soft- reduction, etc., it is necessary to 
give large stiffness, make movable in the longitudinal direction and lockable. Since these mechanisms are available 
and feasible by known technology, the detailed mechanisms are not shown in this description. 

Symbols 124 and 125 are the rolls to give a small amount of reduction gradient. The roll is attached with roll 
crown at the central part to avoid unnecessary and detrimental plastic deformation at the corners of cast piece and 

30 also to effectively transfer the compressive deformation onto the central mushy zone. In the case of this example, 
the reduction is done by oil-hydraulic cylinder 127 that is attached to upper bearing unit. The oil-hydraulic cylinder 
is not necessarily attached at upper side. The plural number of rolls are arranged in the longitudinal direction and 
the reduction force is independently controlled by each roll. The prescribed amount of reduction is to be given by 
the reduction force. Generally speaking, the reduction force needs to be increased as it goes downstream with 

35 thicker solidified layer as shown in Figure 61 . Since the amount of reduction is small (the order of solidification con- 
traction in mushy zone), the stroke of oil-hydraulic cylinder may be small and therefore the length of cylinder may 
be short. Care must be taken on the occasion of a detailed design so that the strokes of both sides of the cylinder 
become equal or even if a little difference occurs, there be no obstacle in the operation. 

Furthermore, the roll is equipped with the driving unit (usually at the lower roll). The number of driving rolls may be 

40 determined by the magnitude of required driving force, etc. The detail of the driving unit is not shown here because 
it can be easily assembled by known technology. 

Next, the relation of the distance between superconductive coils and the width of the coil is stated. It is under- 
stood that the distance between the coils b may be reduced to obtain a stronger magnetic field from above-men- 
tioned Eq. (77). Also, the coils with the relation of a=b is called Helmholtz-type which is possible to obtain highly 

45 uniform magnetic field. 

Specific Example 2: Application to the case where the distance between coils is shortened. 
Considering this point, the width of the coil was expanded and the distance between the coils was shortened to 
obtain a stronger magnetic field compared to Specific example 1 . The cross -sectional view is shown in Figure 63 
(a). For this, the distance between the coils was shortened by setting up spaces in the upper frame 117 and the 

so lower frame 118 to house the roll units. If the distance between the coils is too short, the coils collide with or 
approach too close to the cast piece at the position where the coils cross the cast piece. In such a case, horse- 
saddle type coil may be used to secure a necessary space at both ends of the coil as shown in Figure 63 (b). 

This example is basically applied to the case where the driving units of rolls are not required by properly adjust- 
ing the balance of the Lorentz force and the drawing resistant force due to soft-reduction (Function n of the above- 

55 mentioned). [If the roll drive is inevitably needed, it will become possible to chain-drive toward longitudinal direction 
of the cast piece by attaching gears to the roll edges.] Other mechanisms are similar with those of the specific 
example 1 (omitted). 
Specific Example 3: Application to slab. 
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The specific example applied to a wide slab is shown in Figure 64, where both Lorentz force and small amount 
of soft-reduction are given. Because the reduction rolls are long and slender and easy to bend due to reduction load 
and thermal stress, split rolls are used. The reduction force is given by oil-hydraulic cylinder. The reduction is done 
by upper roll in general and the hydraulic cylinder 127 is attached to each bearing unit. The cylinder stroke may be 

5 short as already mentioned. If more margin is necessary, the unit of specific example 2 can be used. As to the roll, 

one piece of roll may be used whose diameter is squeezed at bearing units or may be divided into split rolls and 
each of them supported independently at the bearing units. Also, it is desirable to prevent the plastic deformation 
at the cross-section corners by attaching roundness at both sides of the end rolls. Roll-drive is done by lower rolls 
in general. Other mechanisms such as electrode are same as those of specific example 1 ; therefore, omitted. 

10 Specific Example 4: Application to parallel casting of plural number of cast pieces. There are two types in the con- 
tinuous caster that cast plural number of cast pieces in parallel at the same time: One is the case that the distance 
between the cast pieces is sufficiently wide and the other is the case of narrow distance. In the former case, the 
electromagnetic booster and the reduction unit may be installed independently. This example refers to the latter 
type. In this case, the neighboring cast pieces are connected by flexible bus bar or cable 131. The electrode box 

75 105 is fixed to the electrode frame 1 07 that extends toward the longitudinal direction of cast piece. Also, plane mill- 
ing tool 108 and gas shield box 109 are attached onto the surfaces of the cast pieces. The magnetic field is gener- 
ated by a pair of upper and lower superconductive coils. The roll reduction units are built in to each cast piece. 
Other mechanisms are as already mentioned. 

20 The case that only electromagnetic force is applied: In this case, the roll reduction function shown in the above four 
examples is unnecessary. Hence, the specific example is not shown. In the case of blooms, etc. mentioned in the spe- 
cific examples 1 , 2 and 4, the cast piece is sustained firmly by the solidified layer. Therefore, upper rolls to support the 
cast piece are not necessarily needed or only minimum number of rolls are required. On the other hand, an appropriate 
number of lower rolls are necessary to support the cast piece, considering that fairly strong electromagnetic force is 

25 exerted onto the cast piece. 

[0180] Considering that, in the case of the slab of the above specific example 3. the width of mushy zone is wider than 
that of bloom and fairly strong Lorentz force acts, it is necessary to firmly support the cast piece by upper and lower rolls 
as in the case of conventional slab casting. 

[0181 ] Also, in the case that the electromagnetic force is strong enough to give rise to an exccessive tensile force in 
30 the solidifying shell with mushy zone both for bloom and billet, a means is possible to relax the tensile force by installing 
a conventional roll reduction unit (not shown) at the downstreamside of the booster of Figure 1 and thereby applying a 
brake by the friction force between the rolls and the cast piece due to the roll reduction. 

Outline of the design of electromagnetic force: 

35 

[0182] Its outline is stated below. The magnitude of the electromagnetic force is set 20 times of the gravitational force 
and the shape of superconductive coil is assumed that of circular for the convenience of calculation (according to Eq. 
(77) and Figure 39 (a)). Take J (A/m 2 ) as DC current density to mushy zone, B (Tesla) as DC magnetic flux density, p = 
7.0 (g/cm 3 ) the density of liquid steel and g r = 980.665 (cm/s 2 ) the acceleration due to gravity, then gravitational mag- 
40 nif ication number G is given by the next equation. 

q . JB(N/m 3 ) 0.1 JB(dyn/cm 3 ) {QQ) 
pg r (dyn/cm 3 ) pg r (dyn/cm 3 ) 

45 

Now, the current density is set J = 5 xl 0 5 (A /m 2 ). The corresponding magnetic flux density becomes B = 2.75 (T) from 
the above equation. 

[01 83] In the above specific example 1 of bloom, take the dimensions of cross-section as 300mm x 400mm, the radius 
so of the coil as a = 0.34m and the distance between coils as b = 0.92m. Then, the required current in the coil becomes I 
= 35431 12 (A) from Eq. (77). Now, when the electric current in a superconductive wire (take square) is set at 2000 (A), 
corresponding number of turns of coil N becomes 35431 12 / 2000 = 1772. Take the cross-section area of a supercon- 
ductive wire (square) as 10 (mm 2 ) (accordingly the current density becomes 200 A/mm 2 ), the cross-section area of the 
coil becomes S = 1772 x 10 (mm 2 ) = 177.2 cm 2 . 
55 [0184] Next, when conducting similar calculation as to the specific example 2 where the radius of coil enlarged to a = 
0.48m and the distance between the coils reduced to b = 0.66m for the same bloom as in the example 1 , the design 
parameters become as follows: 
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Required coil current I = 1,877,224 (A) 
Number of turns N = 936 

(Current density of a superconductive wire = 200 A/mm 2 ) 
Coil cross-section area S = 93.9 (cm 2 ) 

5 

[0185] With regard to the specific example 3 of slab, take the cross-section dimensions same as those of above best 
modes No.2 and No.3 for carrying out the invention (i.e., 220mm thick x 1500mm wide) and take a = b = 0.94m (Helm- 
hoftz-type). Then, the design parameters become as follows: 

io Required coil current I = 2,874,853 (A) 

Number of turns N = 1 ,438 

(Current density of a superconductive wire = 200 A/mm 2 ) 
Coil cross-section area S = 143.7 (cm 2 ) 

15 [0186] The design parameters can be obtained with respect to the specific example 4 similarly, but omitted. The 
design parameters of the above superconductive coils satisfactorily enter into the practical range of use of the NbTi 
based superconductive coil at present; thus, there is no technical problem. Furthermore, it is possible to obtain even 
bigger magnetic field. For more details, refer to a standard text book, for example, Daily Industrial Newspaper of Japan, 
"The Application of Superconductivity", by Hiroyasu Hagiwara. In the above calculations the shape of coil was taken as 

20 that of circular for brevity. Also, the number of pairs of upper and lower coils was taken one, but it is possible to install 
plural number of pairs of coils in order to optimize the uniformity and the strength of magnetic field. On the occasion of 
practical design, numerical analysis of static magnetic field may be conducted by finite element analysis, etc. in accord- 
ance with real configuration of coils by taking these points into consideration. 
[0187] The case of applying the Lorentz force toward opposite direction: 

25 [0188] As aforementioned, these inventors have pointed out that if the liquid pressure drop due to the interdendritic 
liquid flow primarily induced by the solidification contraction within elongated mushy zone in casting direction, reaches 
to the criterion of porosity formation (Eq. (65)), the microporosity forms in between dendrite crystals, and triggered by 
this, high solute concentration liquid in the vicinity of the porosity flows in along V type porosity; thus, resulting in V seg- 
regation bands. The porosity lines up in V character pattern as shown in Figure 12 (b) and the flow occurs along V 

30 bands toward casting direction. 

[0189] Accordingly, the V segregation can be lessened by exerting the electromagnetic force toward the direction to 
impede this flow or the opposite direction with the casting direction. This has been shown in the casting experiments of 
steel bloom (Ref . (9)), utilizing linear motor type electromagnetic apparatus (non-contact type) with no DC current elec- 
trodes attached. 

35 [0190] The electromagnetic boosters described in the above specific examples can also be applied for this purpose. 
That is that the Lorentz force may be exerted toward upstream by reversing either the current direction or the magnetic 
field direction. The Lorentz force is exerted in the range from the upstream vicinity of the position reaching the criterion 
Eq. (65) to the crater end. The magnitude of the reversed Lorentz force may be extremely smaller than 20G of above- 
mentioned computational example, because the purpose is to impede the aforementioned liquid flow. If the electromag- 

40 netic force is too small, there is no impeding effect. If it is too big on the contrary, the high solute concentration liquid 
flows toward the opposite direction and results in inverse V segregation; thus, there is no meaning. The magnitude of 
appropriate electromagnetic force can be easily known by experiment on the real machine. Also, the soft-reduction gra- 
dient may be given in addition. 

[0191] It has been stated in the above Ref. (9) that it is very difficult to apply their inverse method to wide slabs 
45 (because of using linear motor type apparatus). On the contrary, since the electromagnetic booster by this invention 
uses the DC current and DC magnetic field, highly uniform electromagnetic force distribution can be obtained for wide 
slabs as well as blooms (some ingenuity is of course done on designing). Thus, it will become possible to effectively 
impede the liquid flow causing V segregation: However, the microporosity would remain to some extent in this method. 
[0192] In the following, some points about the design for electromagnetic booster will be mentioned including the 
so items so far not referred to. 

(1) Appropriate material may be selected for the electrode such as graphite, ZrB 2 and so on in consideration of 
electric conductivity, wear resistance, etc. 

(2) In the case of the above specific example, taking the contact area of an electrode as 1 00mm x 1 20mm. the cur- 
55 rent of 6000 (A) flows into the electrode and then to the bus bar. Copper plate is usually used for bus bar. The cross- 
sectional area of the plate is determined by using the current density of 3 to 4 (A/mm 2 ) and about 10 (A/mm 2 ) when 
water-cooled. In the above specific examples, L-type and plate-type bus bars were shown. These are not indispen- 
sable matters, but the cable knit with copper wire (having flexibility), etc. may be used for example. These are usual 
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conventional technology, and needless to say, may be devised on the occasion of detailed design. 

(3) The electrode and bus bar need to be fixed firmly because the electromagnetic force is generated onto these 
parts, by the interaction of the current that flows in these parts and the magnetic field. 

(4) Because a large magnetic field is generated in the periphery of the electromagnetic booster, the frames, pillars, 
5 roll units and so on that exist in this space are basically made of nonmagnetic materials such as stainless steel, etc. 

Yet, magnetic material (usually, iron) maybe arranged properly to obtain uniform magnetic field. Besides, the prob- 
lems of the effect of magnetic field on various measuring units and the necessity of magnetic shield can be solved 
by known technology; hence, omitted in this description. 

(5) The upper frame 1 1 7 and the lower frame 1 1 8 in Figures 58, 63, 64 and 65 are not necessarily built as one unit. 
io They may be separately made dividing into the rigid frame to house the superconductive coil, the rigid frame to sup- 
port the roll unit, etc. 

(6) The superconductive wire is generally made of composite material where very fine superconductive wires of 
NbTi and as such are mounted in the copper matrix. The coil is made by winding the wires around bobbin (guide 
jig). The superconductive coil is usually used without iron cores. At present, since the coil needs to be cooled to an 

is extremely low temperature (liquid helium temperature, 4.2K) to hold superconductive condition, the inside of the 
cooling chamber 121 is composed of the combined layers of the liquid helium, vacuum heat insulating layer, liquid 
nitrogen layer, etc. Superconductive technology has already been commercialized to many usage such as particle 
accelerator, MRl, etc. It is expected that once high temperature superconductive material(s) is developed and com- 
mercialized in the future, it will diffuse rapidly in many applications. 

20 (7) It is necessary to know the roll reduction force distribution in order to exert compressive force onto the cast piece 
by means of oil-hydraulic cylinder, etc. and thereby to attach the prescribed soft-reduction gradient (refer to Figure 
61). For this, it is possible to know an appropriate reduction condition via minimum amount of experiments by uti- 
lizing analytical method such as FEM. The analytical method is especially useful in the case that the split rolls are 
used for slab or in the case that the crown rolls are used for bloom. 

25 (8) The air gap 1 16 between the gas shield box and the cast piece is made as small as possible and other parts 
are devised to keep sealed from the atmosphere. One idea- is to arrange something like fine stainless steel scrub- 
bing brush in soft contact with the cast piece at the air gap 116. This enables to save the gas consumption, to hold 
the internal gas pressure in the box to slightly positive and thereby to effectively prevent the re-oxidization. 

(9) As a means to remove the surface oxidized layer of cast piece, various methods are available besides the plane 
30 milling tool such as that of fixing a cutting tool in relative to the movement of the cast piece. 

(10) Atmospheric temperature around the cast piece is raised due to the radiation, heat conduction, etc. from the 
surface of the cast piece. Hence, an appropriate cooling measure such as water-cooling needs to be taken to cool 
the roll bearing unit, oil-hydraulic cylinder, bus bar, etc. 

35 [0193] In the above, the mechanisms utilizing plural number of sliding electrodes and the superconductive coils were 
described regarding the specific apparatus of this invention to exert the electromagnetic force toward the casting direc- 
tion. To say more specifically, by adopting the racetrack type or horse saddle type superconductive coils in accordance 
to the shape of the cross-section of the cast piece, by reducing the distance between the coils as much as possible and 
by adjusting the balance of the distance between the coils and the width of the coils, a highly uniform and strong mag- 

40 netic field can be obtained in the wide space including the cast piece, the rolls, the electrodes, etc. (In this point, it is 
very difficult from the view point of the mechanism to produce a highly uniform electromagnetic force in the mushy zone 
of slab such as shown in Figure 64 or bloom such as shown in Figure 58 by a non-contact linear motor-type electro- 
magnetic generator.) 

[0194] Next, as a means to supply DC current, the method of using plural number of sliding electrodes was described 
45 in this invention. This gives birth to such an effect that enables to optionally control the current density distribution and 
thereby to optionally control the electromagnetic force distribution in the longitudinal direction of cast piece. 
[0195] As to the case that the soft-reduction is used as a supplemental means to reduce the required electromagnetic 
force, it is possible to optionally control the reduction force distribution by introducing the independently controlled oil- 
hydraulic system of this invention and thereby to control the gradient of the reduction quantity. In this case, each hydrau- 
50 lie system needs not to be controlled independently for each roll, but may be controlled every several rolls depending 
on the case. Also, oil may not necessarily be used as a pressure transmission medium. Furthermore, by installing driv- 
ing unit to rolls, it becomes possible to adequately control the tensile force that takes place in solidifying shell to prevent 
cracking. 

[0196] As above mentioned, the apparatus by this invention enables to obtain desired electromagnetic force, magni- 
55 tude of reduction force and reduction distribution. As a result, it can realize to obtain the cast piece without the internal 
defects, and also enhance the productivity by high-speed casting. On the occasion of the application, all the functions 
mentioned in this description regarding the electromagnetic force, soft-reduction gradient, control of roll-drive, etc. may 
not necessarily be used. 
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Application Range of the Electromagnetic Continuous Casting Method 

[0197] By the above four best modes for carrying out the invention, the electromagnetic continuous casting method 
by this invention (called E process) was verified, and the specific examples of the electromagnetic apparatus were 

5 shown. The E process can be applied to all continuous casting processes besides the vertical -round type bloom, the 
vertical-bending type slab and the bending type bloom taken up in this description: Namely, recently notified thin slab 
casting (thickness of the slab is the order of 50mm or 60mm at present), so-called near-net-shape casting with irregular 
cross-sections like H-type, and further the composite casting process where different grade of steels are cast simulta- 
neously into the same mold, in addition to the conventional processes such as vertical-bending blooms and billets and 

io bending slabs and billets. The reason for this is that the principles of E process, that pays a special attention to the pres- 
sure drop toward casting direction in interdendritic liquid within the mushy zone at the cross-section of cast piece, holds 
the liquid pressure higher than that of the critical pressure of porosity formation and thereby enables it possible to make 
castings with essentially no central defects (microporosity and segregation), possess an universality to all these proc- 
esses. 

is [01 98] The interdendritic liquid flow toward casting direction induced by the solidification contraction is a physical phe- 
nomenon generally common to metallic alloy. Therefore, E process can be applied to all steel grades such as carbon 
steel, low alloy steel, stainless steel, etc: E process can also be applied to non ferrous alloys such as aluminium alloys, 
copper alloys, etc. 

[0199] E process consists of the method of exerting Lorentz force alone and the combined method of Lorentz force 
20 and soft-reduction. In either case, there is no effect if the timing of solidification, i.e., the position of Lorentz force appli- 
cation (the distance from meniscus) is mistaken. For example, if the Lorentz force is applied at the downstream side of 
the position where the liquid pressure reaches critical pressure of porosity formation, it acts to further accelerate the 
liquid flow to form V segregation because the V pattern of porosity has already been formed. Therefore, there is a pos- 
sibility to create conversely severer V segregation depending on the magnitude of the Lorentz force. If the applied posi- 
25 tion is too upstream side on the contrary, the Lorentz force acts to uselessly increase the liquid pressure where pressure 
rise is not expected. And the effect to the crater end vicinity where liquid feeding is most needed becomes small: thus, 
unfavorable. Even in the case that the applied position is appropriate, if the Lorentz force is too small and become less 
than the critical pressure, there is a possibility to accelerate V segregation due to the formation of porosity. Accordingly, 
it is very important to quantitatively know the position of the critical liquid pressure and required Lorentz force: For this, 
30 the numerical analysis by computer disclosed in this description provides with the most effective means. (It would prob- 
ably be impossible to know this critical position directly by physical measurement. Further, it is impossible to know the 
required Lorentz force distribution by experimental means.) This is the reason why the above-mentioned numerical 
method that these inventors developed is indispensable as a means to comprise E process. 

[0200] This computer program takes the formality to be stored and sold in the forms of sauce program / application 
35 program to the various memory media such as MT (magnetic tape), floppy disk, CD-ROM, DVD, semiconductor mem- 
ory cards, media on the internet, etc. Thus, it is possible to perform a series of analytical work of the input of operating 
conditions, the implementation of computation, display of the results on the computers such as personal computer, 
work station, large frame computer, supercomputer, etc. 

40 Drawing of Cast Piece by the use of Electromagnetic Force 

[0201] The Lorentz force generated in E process can be utilized as a drawing force of cast piece. In bending-type or 
vertical-bending-type continuous casting, the cast piece undergoes drawing resistance such as the drawing resistance 
due to unbending, frictional resistance between the cast piece and the mold, etc. For example, it has been reported in 

45 Ref. (3 1 ) that 60 tons of drawing resistant force was measured on a real machine for slab casting (cross-section 1 90mm 
x 1490mm, casting speed 1.5 m/min). As the casting speed is increased, the drawing resistant force is increased. In 
order to obtain the drawing force corresponding to such a large drawing resistant force, it is necessary to effectively act 
the driving torque by rolls on the cast piece. Hence, a mutti -drive-system is adopted in general. However, it is consid- 
ered that some influences occur to the quality of the product in the method of applying the frictional force onto the cast 

so piece: For example, if the compressive force is too big, the solidifying shell deforms resulting in internal crack and seg- 
regation (Ref. (31)). 

[0202] On the other hand, the Lorentz force acts statically on the cast piece, enables to reduce the number of driving 
rolls and thereby to reduce the above mentioned compressive force. Hence, it will work favorably on the quality. 

55 Utilization Procedure of the Electromagnetic Continuous Casting (E Process) 

[0203] Utilization procedure of E process to real continuous casting is as follows. 
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(1) Matching the numerical analysis with the test on a real continuous caster. 

[0204] The computational results described in the above best modes for carrying out the invention are of course 
accompanied by errors. The primary cause of the errors is associated with the heat transfer coefficient on the surface 

5 of cast piece and the accuracy of various physical properties. With respect to the physical properties used in this 
description the values are reasonable that are quoted from various published references, but it is difficult to expect 
accuracy for many data. The second cause of the errors lies in the modeling of the morphology of dendrite crystal and 
the accuracy of resulting permeability K. The validity of the modeling of complicated dendrite morphology has been 
proved by Ref. (18). In addition, it is known that the permeability K differs in the parallel direction to the growth direction 

io of dendrite crystal (Kp) and in the perpendicular direction to it (Kv) (Ref. (32)). It seems that Kp and Kv depend on the 
cooling rate. However, the reality is that there is no accurate data on the relationship between Kp and Kv of commercial 
steels. Therefore, upon the matching, above two points needs to be taken into account. 

[0205] It is possible to correct the error by the above primary cause by measuring the surface temperature (or internal 
temperature) change (for example, Ref. (33)). Good amount of data have been accumulated so far about the relation- 
15 ship between the cooling conditions such as water-mist spraying and the surface heat transfer coefficient. Now that the 
measurements o1 solidifying shell thickness and crater end position have become possible, an accurate correction is 
possible. The method of correction is optional. To cite an example, the correction can be done by thermal diffusion coef- 
ficient X7cp. 

[0206] With regard to the error by the second cause, the error may be corrected so that the calculated and the meas- 
20 ured values of the critical position of porosity formation coincide, introducing a parameter a K to correct the influence of 
anisotropy of columnar dendrite into the equation of permeability K (Eq. (27)) besides the correction factor a in Sb equa- 
tion (Eq. (28)). That is, the critical position of the porosity formation may be designated by observing the conditions of 
internal defects (the range of formation, the size of porosity, etc.) and comparing with the numerical results. 

25 (2) Determination of the optimized condition of E process via numerical analyses 

[0207] Once these correction factors are established by the above procedure (1), the optimized conditions to eliminate 
the internal defects can be found by fully utilizing the numerical computations: i.e., the range and the magnitude of 
Lorentz force, soft-reduction conditions (if required), etc. This has already been described. 
30 . [0208] Since thus determined optimized conditions are those corrected in the procedure (1), they are highly reliable. 
On the occasion of real operation, these parameters are set to safer side, needless to say. 

[0209] This invention is composed and works as above mentioned, enables it possible to predict the position, the 
quantity and the range of internal defects of continuous castings, and is able to evaluate optimum applied range and 
magnitude of the Lorentz force to suppress the formation of the internal defects. Thus, this invention can provide with 
35 unprecedented, excellent method and apparatus for continuous casting which enables it possible to obtain highly qual- 
ified continuous castings with essentially no segregation or no porosity regardless of the chemical compositions. 
[021 0] Because this invention can combine the electromagnetic force with the soft-reduction, it becomes possible to 
obtain highly qualified steels with essentially no segregation or no porosity in high-speed casting as well. 
[021 1 ] At last, the effect of this invention can be summarized as follows: 

40 

(1) The internal defects (central segregation and porosity) can completely be eliminated. 

(2) High-speed casting is enabled. 

(3) The degree of freedom for chemical compositions can be expanded. 

(4) The variety of steel grades of continuous castings can be expanded. 
45 (5) Drawing apparatus can be made simple and effective. 

[021 2] Especially, with respect to the above (2), the number of continuous casters can be cut to half by increasing the 
casting speed 2 or 3 times. Its economic effect is significant. As the magnetic apparatus, superconductive magnet is 
preferred to conventional electromagnet from the viewpoints of construction cost, operating expense, energy conserva- 
so tion and saving of the space. 

[021 3] thus, the continuous casting process by this invention can be said a new process excellent in productivity and 
economic efficiency as well as in the quality. 

[0214] These inventors, without staying at the macroscopic phenomena such as heat and fluid flow, have coupled 
these macroscopic phenomena with the microscopic solidification phenomena such as dendrite growth and solute 
55 redistribution in multi-alloy system, and developed the computer program where the effects of electromagnetic force, 
mechanical deformation and pearlitic transformation were further incorporated. To the best of these inventors' knowl- 
edge, the whole picture of the internal defects problem in continuous casting has been grasped for the first time. 
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Discretization of the Governing Equations 
A: Discretization of the Energy Equation 
s [0215] The discretization equation regarding temperature is as follows. 

a p T p = a N T N + a s T s +a T T T + a B T B + a w T w + a E T E + b 

(A.l) 

10 

a N =[D„AQP n \) + {-F„,0)]A„ (A .2) 
a s =[D s A{\P s \) + (F s ,0)]A s (A .3) 
^KA^D + H^O)]^, ( A. 4) 

o fl =[D^(|/i|) + (^,0>]^ (A. 5) 

= t A„ ^ (I y w | ) + (- f v „. ,o))a w ( a. 6 ) 



25 



30 



35 



40 



45 



50 



55 



a E =[D e A(\P e \) + {F e ,0)]A e (A . 7) 

a / J - fl /V + a 5 + a r + ^5 + a W + a E + a °P ( A. 8 ) 

old f ~~\old AK 

=\CP)P — (A..9) 
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10 



b = a?r p u + \L(p)f( gs - gf)„ + - P°'% }^ 

A/ 

^^AK+^L + (4-^)p 5 r}; W [v 5 • Vg s ] 

+ " 4 )p s T% d g SJ ,\y • Vj - (c* A )f &i p[V s • vn 

tv s • %sJ=v£p(gs < n^s / s) A p +v |p<gs,t-gs,b)At+vf / p(gs,w-gs^) A w 



(A. 10) 
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(A. 11) 


[V • 


V S ]=A n vS n - A s vf # s + A t (v| t -v| b )+ A w 


(v |w- v 3,e 








(A. 12 ) 


[v s • 


VT]=vf p (T n -T s )A p +vf p (T t -I^)A t +vf^ p (T w 




(A. 13) 


D n = 


~K< S n i ^ =i C PPLSL V \) n 




(A. 14,15) 


A = 


Z s /S s ;F S =(cpp L g L v i ) s 




(A. 16 , 17 ) 


A = 


J t / 6, \F t ={c ] £p L g L v 2 ) [ 




(A, 18,19) 


A = 


A b /5 b ; F b ={CpP L g L v 2 )b 




(A. 20,21) 


A,= 


■ 5 W ; F w = {cpp L g L v 3 ) w 




(A. 22, 23) 


A = 


J e / 5 e ; F e =(cpP L g L v 3 ) e 




(A. 24, 25) 



Where, Peciet number P n = F n /D n , etc. 
Function A(\P\) = <0,(1-0.1 • |P|) 5 >, etc. 

Symbol < > means to choose the bigger of the values in the bracket 
40 The lower suffixes 1 , 2 and 3 of the velocity express the velocity component in N, T and W directions respectively view- 
ing from the grid point P, and are defined at the face n, s, ...of the element (control volume). Upper suffix old denotes 
the value at the time At before. Also, X takes the harmonic mean of the neighboring two elements at the element face. 
Namely, 

K= . 8 " ; .ate. 
8n" is the distance between P and n, and 6n + n and N. 

50 

B: Discretization of the Solute Redistribution Equation 

[0216] The discretization equation of the liquid solute concentration C n (written as C for brevity) is as follows. 

55 
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a p Cp = a N C N + a s C s + a T C x + a B C B + a w C w + a £ C E + b ( B . 1 ) 



10 



15 



20 



25 



30 



35 



a N =[D„A(\P„\) + (-F n ,0)]A n 
a s =[D s A(\P s \) + {F s ,0)}A s 

a r =[D t A{\P t \) + (-F t ,0)]A t 

a B =[D b AQP b \) + (F b ,0)]A b 

aiv=lD„AQPJ) + (-F„,0)]A w 

[D e AQP e \) + (F e ,0)]A e 
a N + a s +a T +a B +a w +a E + a P 

A V 



a E 
a p 



At 



b=4- 



i (^° ,d 



V J *\ ) 

f \oU , / 'old 



a +a k-vd] 



(B.2) 
(B.3) 

(B.4) 

(B.5) 

(B.6) 

(B.7) 
(B.8) 

(B.9) 



(B. 10) 



Where. A n is given by Eqs. (12) and (16), B n by Eqs. (13) and (17), C „ by Eqs. (14) and (18) and D „ by Eqs. (15) 
and (19) in this description. Upper suffix * in the term b expresses the updated value in the iterative convergence com- 
40 putation and takes the average value of the time t and t-At (Crank- Nicholson scheme). 
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[v- (&v s )]=(&vf )„4, -C&vfH -(&v* 

For alloy element n, 

v s ■ vg]=vfp(cf v ,-cf, 5 )4 i3F +^(cf a -^)4 +v^(cf, w -^e)A 

In the above equation, C y is used instead of Vf, for j-type alloy (equilforium solidification). 



(B.ll) 



(B.12) 



48 



BNSDOCID: <EP O90959BA1 J_> 



EP 0 909 598 A1 



15 



D n = D n / 5 n \F n =v l n (B.i3) 
D s = D s i 5 S ; F s = v, s (B. 14 ) 

D { =Djl 5 , ; F, = v 2f , (B.i5,i6) 

-Oft/ $ b> F b = v 2 : b (B. 17,18) 

w=DJ $ n> ; F u . = V 3?w (B.19,20) 

D e = Z) e / £ e ; F e = v 3 . (B.21,22) 



Peclet number P n = F n /D n , etc. 
20 Function = (0,(1-0.1 • |F|) 5 >, etc. 

The meanings of symbol < ) and lower suffixes 1 , 2 and 3 of velocity are the same as in the previous section A. Also, the 
diffusion coefficient D (= D 0 exp (-Q/RT)) similarly takes the harmonic mean at the element face. There are as many 
equations as the number of alloy elements. 

25 C: Discretization of Temperature- Volume Fraction of Solid Equation 

[021 7] The discretization equation of temperature T vs. volume fraction of solid g s is as follows. 
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50 



a P T P = a N T N - + a s T s + a T T T + a B T B +a w T w + a E T E +b ( c.i) 
a N =(-F n ,0)A n ;F„ = v ln (c.2) 
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w 
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45 



a s =(F s ,0)A s ;F s =v u 



a T ={-F t ,0)A, ; F t = v 2 , 
a B=(f r b>°) A b ; F b= v 2,b 
a w={~F w ,0)A w ;F W =v 3w 
a E =(F e ,0)A e -F e 



— v- 

e i.e 



a P ~ a N +a S +a T + a B + a W +a E + a P 



a D = 



AV 
At 



= ld + Sl + S2 | + S3 + S4 



(C.3,4) 
(C.5,6) 
(C.7,8) 

(C.9,10) 
(C.11,12) 

(C.13) 
(C. 14) 

(C.15) 
(C.16) 



C. 17 ) 



c^yT, ^ ,v ni L H^^m 

3 



j] 



iy-(9B v s> ^ 



s 4 =- 



f 



— s 

n 



(c.18) 

(C.19) 



Where [v • (g s v s )] is given by Eq. (B.1 1). Because the influence of S 4 is small, it was neglected. A n is given by Eqs. 
so (12) and (16). B n by Eqs. (13) and (17), C „by Eqs. (14) and (18) and 6 „by Eqs. (15) and (19) in this description. 

D: Discretization of the Darcy Equation- Pressure Equation 

[0218] the velocity equations by Darcy's law (Eq. (26) in this description) are as follows. 
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\n = 



f 



GF Ul+ EMF , „+ — ~ /J,V 



S n 



V 2,< = 



10 



GF 2t +EMF 2t + 



Pr 



St 



V 3,n. = 



75 



KPSlJ. 



ow 



(D.l) 



(D.2) 



(D.3) 



Where, GF-j, etc. and EMF^, etc. are the components in X 1t X 2 and X 3 directions of gravitational and electromagnetic 
forces, respectively. Namely, GF^ = aGF^p L g n etc. 

aGF 1( etc. are the coefficients in X^... directions of the curved coordinates (X 1t X 2 , X 3 ). For example, aGF 1 = aGF 3 = 
20 0 and aGF 2 = -1 for a vertical continuous casting. Suffixes 1 ,2,3 of K take into account the anisotoropy of columnar den- 
drite crystal: For example, in the case of slab, take K 1 = K P (parallel to the growth direction of columnar dendrite) and 
K 2 = K 3 = K v (vertical to the growth direction). In the case of equiaxed crystal, set K-j = K 2 = K 3 . Also, the harmonic 
mean is taken at element's faces. 

[0219] Next, the pressure equation is derived by combining with the continuity equation (Eq. (9) of this description and 
25 the above (D.1)). 

Thus, discretization of Eq. (9) is given by, 

— — old AV , ^ 

(P-P )p — + (PlSl v xXA -(PlSl>\\A +KPLg L v 2 \ -(p L gS' 2 )M 



30 



35 



A(PLgL V ^-(pLgL V 3 lK+(Psgs^)nA-(Psg S ^)sA 



(D.4) 



On the other hand. Eq. (D.1), etc. are put into the following form. 



40 



(JL) 

U&J„ 



(JL) 



IS 



(D.5,6,7) 



45 



SO 



Km). 



{GF^EMF^)-d = 



ISe 

(D. 20,21, 22) 



Substituting these equations into Eq. (D.4) and arranging, the following discretization equation is obtained. 
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apPp = a N % + a s % + a T Dp 4- a B Pfe + a w % + a £ P£ + b 


(D.23) 


5 


a N = (P LSl)n d n A n 


(D.24) 




a S = (P LSL)s d s A s 


(D.25) 


10 


a T = (P LSL)t d t A t 


(D.26) 




a B = (P L§L)b d b A b 


(D.27) 


15 


a \Y= (P LgL\r d u- A M- 


(D. 28) 




°£ = (P L2L)c d eA c 


(D. 29 ) 


20 


cip = citf + a$ + cip + a $ + a w + 


(D.30) 


25 


b = Cp° ld p) p A V + [^vL] + [rfrvS] 


(D. 31) 
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40 



A A 

[divL ]=A s (p L g L vi) s -A n (p L g L vi) 



n 



A A 

(PL9L v 2)b-(PLgL v 2)t 



+ A w | 



A A 
(PLgL v 3) e -(PLgL v 3) w 



(D. 32) 



45 



SO 



[divS] - A s (p s g s vf ) s - A n (p s g s vf) n 

+ A t{ps9S v | ) b " (Psgs v |)t} 
+ A w Jp s g s v|) e - (p S gsv|)w} 



( D . 3 3) 



55 



Where, p = p L g L + p$g s and g L + g s + g v = 1 . It is 1 Semphasized again that P field is determined so as to satisfy the 
continuity condition which include the effect of Lorentz force, in addition to porosity, solid deformation and gravitational 
force. 

E: Discretization of the Momentum Equation 

[0220] Staggered grid is used for the momentum equation (refer to Ref. (20)). The discretization equation of v n is 
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expressed as follows by the use of staggered grid in (r) direction (refer to Figure 41 (a)). The suffix 1 is omitted for 
brevity. 
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a n = a nn + a s + a NT + a NB + a NW + a NE + a » ~ S >- 



VI ff 



^ = [d,a<\p n \) + (-f,.,o)]a, 

a s =[D P A(\P P \)+(F p ,0)]A p 
a yr =[D nl A\p n ]) + (-F nn 0)]A nl 

^ = k*^|)+(^,o>K* 

«v 1( =k,^^,,,|)+(-^^o)K,,. 

^=k^^|K(-^ t oO)K t . 

o , -old AV 

^,=(PLgL+P -P)n— ~ 

A / 

6 = £1^+^+^.^.09^3)1, 



(E-l) 

(E.2) 
(E.3) 
(E.4) 

(E.5) 

(E.6) 
(E.7) 

(E.8) 
(E.9) 

(E. 10) 



(Symbol * of v n denotes the updated value in iterative convergence computation. Similarly thereafter) 

(E.ll) 



Old ( n n \<>ld 



A; 



40 



45 



SO 



As for orthogonal coordinates system (Cartesian coordinates): 



S c _„={GF x +EMFO n -&V n 



(E.12) 
(E. 13) 



[V-( P p S V S )] n ^(p^ S vfW%,(R^ S vf)p+ 

(E.14) 



55 As for the staggered grid in r direction of the orthogonal curvilinear coordinates of Figure 9: All other terms are common 
except for S c n and S n , which are given as follows. 
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(E.15) 



10 



Ax,-(p L v 2 ) n Ax x A0Ax 3 ( E - 16 > 



T5 As for the staggered grid (omitted in this description for want of space) in r direction of (r.O.z) cylindrical coordinates. 



20 



A0Az 



S Cin is same as Eq. (E.13). 

D N =fi N ISN;F N ={p L g L v\) N 
D p =m pi 8 P\F P = {p L g L v { ) P 
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D nt =M nt /<5 NT \F nt ={p L g L v 2 )nt 

Dnb=Vnb ISNB > F nb = (PlSl V 2 )nb 
A»v = M nJ 5 NW i F nw = (PlSl v 3 ) 

D ne = M ne l S NE \F ne = {p L g L v 3 ) n 



(E.17) 



(E. 18,19) 
(£.20,21) 

(E.22, 23) 

(E. 24,25) 
(E. 26,27) 
(E .28,29) 



n takes harmonic mean. 

45 Using staggered grid in 2 direction (refer to Figure 41 (b)), the discretization equation of v 2 is as follows (suffix 2 is omit- 
ted for brevity). 



so 
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10 



is 



so 



25 



30 



a 



a 



a , v t = a nnr V nnt +a ssKsr +a t t v t ( +a b V b +<W v h-my +a eet v eet + 

+{P p -P T ).A t 
a t = a nn t + a sst +a tl +a b + a ww[ + a ee t + a °- S, 

= [D nl A\p„])+(-F„ n 0)]A nl 

a tl = [D T A§P T \)+(-F T ,0)]A t 

a b =[D p A(\P P \)+(F P ,0)]A t 

a ,=[a,,^^,,|)+(-^ 5 0)]^ m , 

^=[A,^|)+(^»0)l^-r 



a»={p L gL+P° ld -p) l 



A / 



6 = a? w v? /rf +1 S c ., +v;[V.(p A v s )] ( 



a t = (PlSl ) / — 



(E . 30 ) 

(E.31) 
(E.32) 
(E . 33) 
(E.34) 
(E.35) 

(E . 36) 
(E.37) 
(E . 38) 

(E.39) 
(E . 40) 
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[V • (p S g S v S>] t = A nt (p s g s vS ) nt _ AstCpsgsvf ) st 

+A t {[p s g s v|) T - (p s gsv|)p} 

+ A w4psgs v |)wt - (PsgsvfJet} 



(E.41) 



As for Z direction staggered grid of orthogonal and cylindrical coordinates: 



50 



55 



55 



BNSDOCID: <EP 090959BA1_I_> 



EP 0 909 598 A1 



s=- 



V K 2 J, 



AV, 



Sc., = (GF 2 + EMF 2 ) r AV t 



(E.42) 



(E.43) 
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As for Z direction staggered grid of curvilinear coordinates system (Figure 9): all common except that S c t and S t are 
given by the following equations. 



l K 2j t 



Av t~ (pLgL)t(- v i, t ,oW 1 Ax 3 A0-^ t ln 



v r nty 

(E.44) 



/t =(GIfe+EM%> t . AV t +(p L g L ) t (| v * t ,o) V ^ t Ax 1 Ax 3 AG 



2 Mt< v l,T- v l,p) ln 



Ax 3 



A,/ Snnt;F nt = (p^v, ) M/ 

A/ =M st l 5 sst\F st =(p L g L v l ) sl 
D T = i u r /S tt;F T = (p L g L v 2 ) r 
Dy> =M P I S b\F P =(p L g L v 2 )p 

D wr =M wl !SW;F w ,=(p L g L v 3 ) wt 
D et = fi et l 5E\F et =(p L g L v 3 ) et 



(E.45) 
(E. 46,47) 

(E. 48,49) 
(E.50,51) 
(E. 52,53) 



(E .54,55) 
(E. 56,57) 



so Using the staggered grid in X 3 (Y) direction (refer to Figure 41 (c)), The discretization equation of v 3 is as follows (suffix 
3 is omitted for brevity). 
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a w = a NW + a SW + «7W + a BW + Oh W + ^ + «° ~ 5^. 

^;'=h^j)+(-^,0)k'- 
^=k,^,|)+(^,0)}4, M , 

^,'=k^|^,|)+(-^ 3 0)}4 ni . . 
a ww = |pw A (Nto|) + (- f w,o)]a w 

= [d p a(\p f \) + (f p ,o)]a w 

o , -otf — AK 

a w={PLSL+P -P), 



A / 



* = + 5 e . lr + v;[v.(^ s v 8 )i r 



A/ 



+^v|(P^S v 2)tw-(P^S v 2)b W } 



(E.59) 

(E. 60) 

(E . 61 ) 

(E . 62) 

(E. 63) 
(E . 64 ) 
(E.65) 

(E.66) 

(E.67) 
(E .68) 



(E.69) 



As for Y direction staggered grid of curvilinear (Figure 9) and orthogonal coordinates systems: 



s.= 



V K 3 J 



AK.. 



S c ^ =(GF 3 + EMF 3 ) W AV W 



(E.70) 
(E.71) 



57 



BNSDOCID: <EP 0909S9BA1_I_> 



EP 0 909 598 A1 

As for e direction staggered grid (omitted for want of space) of cylindrical coordinates (r. 0, z): 
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s = 



K K 5 J 



(E.72) 



S CtW = (GF 3 +EMF) W AV W +(p^J^(-v i ;.,0)v: u .^A^ 



+2 ^u.(v UI -v lrP )ln 



f r \ 



■Az 



D nw = P n J S NW ■ F n w = (p L g L v x )„ ». 

A, = P S JSSIV;F SW ={p L g L v l ) sw 

A =V,J5 TW ; F t w = (p L g L v 2 ) , w 
F> bw = P bw I 5 BW ; F bw = (p L g L v 2 ) bw 

D w = M wl 5 \vw ; F w = (p L g L v 3 ) w 
D P = fi p l Se ;F P = (p L g L v 3 ) p 

Pressure discretization equation: 

[0221] The discretization equation regarding the pressure in the momentum equations Eqs. (E.1), (E.30) and (E.58) 
are derived by combining these momentum equations with the continuity equation (Eq.(9)), similarly as in the case of 
Darcy analysis. For this, the momentum equations are deformed as follows. 



(E.73) 

(E. 74,75) 

(E. 76,77) 

(E. 78,79) 
(*E. 80,81) 

(E. 82,83) 
(E. 84, 85) 
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= TfljjPy* *± + P £ -P i +ds( r s - Pp) (E . 87) 

= Z<W^+* + Pp^L . Ai = v 2> , + d,(P P -P T ) (E . 88, 



i(=£! ^ +!m , 4 _ ^ ^ _^ (E . g9) 





T ■ 

a, 



























4 V = V 3 .w +d w (P p -P w ) ( E . 9 0 ) 



^vv=\e + <(^-^) (E.91) 



[0222] Symbol £ is the sum of the products of the coefficients and the velocities of surroundings. Substituting the 
above Eqs. (E.86) to (E.91) into Eq. (D.4) and arranging with respect to P, the discretization equation of P is obtained. 
P is, as seen from Eq. (E.92), defined not in the staggered grid but in the original grid (refer to Figures 8 and 9). 

a P P P =a N P N +a s P s +a T P T + a B P B +a w P fV + a E P E +b (E.92) 
a F = a N + o s + a f + a B + a w + a E ( E . 9 3 ) 

aN=<J>L8L)n d *A n '> d n = *I>N /a n (E. 94,95) 

a S ={PLS L )s ci ^s \ d s = A ps /a s (E. 96,97 ) 

ar=iPLSL)td,A t \ d t =A t /a t (E.98,99 ) 

a B =iPLSL)b d b^t \ d b = A t /a b ( e.ioo,ioi) 
air = (p L gL) w d w A w ; d w = A w I a w <e.102,103) 
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a E = (PLgL)e d e A u -i d g = A w / a e (E. 104, 105) 

— old — £j' r 

b=(p -p) + ( PLgL v i ) A -ip L g L \) n A +m&yz) b -(p L gA),H 

01 (E.106) 
+{(^&vg e -( Pt& v 3 ) > X-[V(p^v s )] 



10 



The term - [ v - (ps0 s v s )] in Eq. (E.106) is given by Eq. (D.33) (note the negative sign in front of bracket [ ]). Pressure 
correction and velocity correction equations: If the velocity field is converged in iterative computations, the pressure field 
can be obtained at once by solving the pressure equation. Thus, it is required to correct the velocity field. For this, cor- 
15 rect the pressure field first. This is the iterative solution method by SIMPLER algorithm. The algorithm is as follows: now, 
put as follows. 

P (solution) = P* (updated value) + P '(correction) (E.107) 
20 V (solution) = v* (updated value) + v f (correction) (E.1 08) 

The momentum equations for P and P* are given respectively as follows. 

°n v Ln = T a nb v Lnb + b + (Pp ~P N )A pN (E.109) 
a n v l,n =X a nb v \nb + b + ( P P ~ P n) A PN ( E . 1 1 0 ) 

30 

Take the difference between the above two equations and regard as follows (for the sake of convenience). 

^ n ={Pp-P N )— = d n {p P -P N ) (E.lll) 
40 a n 

By substituting Eq. (E.1 11) back to Eq. (E.108), a series of velocity correction equations are obtained as follows: 

45 



50 
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( E . 114 ) 



( E . 112) 



(E . 113) 



(E . 115) 



(E . 116 ) 



>3 f e =V 3 e+d e (P E -P p ) 



(E.117) 



d n . • • • are given by Eq. (E.95) • • Substituting Eqs. (E.1 12) to (E.1 17) into the continuity Eq. (D.4) and arranging 
with respect to P\ the pressure correction equation is derived as follows. 

a P Pp = a N f£, + a s P' s + a T F^ + 3l B Pq + a w P^ + a E Pfc + b (E. 118) 

Where the coefficients a P and a N , • • • are given by Eq. (E.93) and Eq. (E.94), • • • , respectively, b is given by Eq. 
(E.106). Yet, v\ n , - • • are used instead of v 1r) , • • • . Also, when v* field is converged, b = 0; therefore, the con- 
vergency is judged by if b« 0 (a small number). 
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Table 1 Explan ation of symbols used in this 

description 

Energy equation: 



T Temperature (°C) 

t Time ( s ) 

io At Time increment in computation ( s ) 

L S o 

Cp^Cp Specific heat of liquid and solid (cal/g°C) 

Pl>P$ Density of liquid and solid (g/cm3) 

75 ^L^S Thermal conductivity of liquid and solid 

(cal/cm s°C ) 

g s Volume fraction of solid 
2Q g L Volume fraction of liquid 

gy Volume fraction of porosity 

p Average density (g/cm3) of mushy zone: 

25 given by PsSs+PlSl 
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A Average thermal conductivity (cal/cms c C) 

of mushy zone : 

given by ^sSs+ ^ lSl 
V/, Liquid flow velocity vector (cm/s) 

V S Deformation velocity vector of solid 

(cm/s ) 

L Latent heat of fusion (cal/g) 

Qj Heat of Joule by electric current 

(cal/cm3s) 



Solute redistribution equation : 
_ 

Liquid concentration of solute element n 

(wt%) 

Average solid concentration of solute 
element n (wt%) 
Cn Average concentration of solute element n 

in mushy phase (wt %) 

Diffusion coefficient of solute element n in 
liquid ( cm^/s ) : 

given by D L = D Q exp(—Q/ R T) 

D 0 = Diffusion constant (cm 2 / s ) 
Q = Activation energy of diffusion 
( cal/mol ) 

R = Gas constant, 1.987 (cal/molK) 
T = Absolute temperature (K) 
C n Solid concentration of solute n at solid- 

liquid interface (wt%) 

L S 

m n r> m n Slope of liquidus and solidus line of 

solute element n (°C /w t %) 
T. x( Melting point of base metal (°C) 

CrO 
^ Solute content of alloy element n (w t%) 

T n Starting temperature of solidification for 

C n in an alloy of base metal-alloy element n 

phase diagram. 

old 

Ss Volume fraction of solid at time t-A t ( s ) 

rr° ld 

Sv Volume fraction of porosity at time t-A t 

( s ) 
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Solidification contraction defined by 

{Ps-Pl)'Ps 



Darcy equation and Momentum equation 



JJ Viscosity of liquid (dyn s /cm 2 ) 

K Permeability (cm 2 ) 

P Liquid Pressure (dyn/cm 2 ) 

X Body force vector (dyn/cm^) 

g r Acceleration due to gravity 980. 67 (cm/ s 2 ) 

Sb Specific surface area of dendrite crystal 

(cm 2 /cm3) 

f Dimensionless constant in permeability K 

having the value of 5.0 
<j) Configuration coefficient of dendrite cell. 

2/3 for cylinder, 
d Diameter of dendrite cell (cm) 

° L.S Surface energy at solid-liquid interface 

(cal/cm 2 ) 

Electromagnetic analygj? 

f Electromagnetic body force vector (Lorentz 

force) (N/m3) 

J , J Current density, Current density vector 

(A/m 2 ) 

B Magnetic flux density vector (Tesla) 

C Electric conductivity (1/Qm) 

E Electric field vector (V/m) 

<j) Electric potential (V) 



Table 2 Physical properties of IC-lCr steel and 0.55% 
carbon gteel 

IC-lCr bearing steel 0.55%C steel 

Specific heat of liquid C L (cal/g°C) 

0.15 0.158 
Specific heat of solid Cg (cal/g°C) 

0 . 15 Figure 25 
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15 



Thermal conductivity of 1 iquid X L ( cal /cms c C) 

0.083 0.071 
Thermal conductivity of solid X s (cal/cms c C) 

0.064 Figure 25 

Solid density ( austenite ) p s (g/cm 3 ) 
10 7.34 7.30 

Density of pearlite p p (g/cm 3 ) 

7.8 7.8 
Latent heat of fusion L (cal/g) 

66.0 65.0 
Viscosity of liquid ^ (poise) 

0.085 0.08 
20 Radiation coefficient at surface £ 

0.3 0.3 
Latent heat of pearlitic transformation Lp (cal/g) 

20.0 20.0 
Upper temp, of pearlitic transformation (°C) 

735.0 760.0 
Lower temp, of pearlitic transformation °C) 
30 400.0 300.0 

Surface tension at liquid-CO gas interface 

a LV (dyn/cm) 1700.0 1700.0 



25 



35 



40 



45 



Constant in liquid density by p L = p \+ 7] h n C„ + h° 



L 



Constant p L ( g / cm 3 ) = 9.265 
Constant h° ( g / cm 3c C) =-1.45xl0**3 

h n (g /cm 3 • wt%) Equilibrium partition ratio 
C —0.08 Locally linearized (Figure 14) 

Si —0.087 0.5 



Mn —0.014 0.75 

Cr —0.059 0.85 

Ni 0.004 0.95 

50 P -0.084 0.06 

S -0.09 0.05 
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Physical properties in Sb equation Eq.(28) 
Configuration factor of dendrite (j) = 0.67 
5 Solid-liquid interfacial energy o LS ( C al/cm 2 ) = 6 x 10~ 6 

D q and O in diffusion coefficient 
D L ~D£exp(-Q/ RT) 

■ — ^ 0 l crn L s > activation energy Q (cal/mol) 

C 1.74 x 10-3 7570 

Si 7 . 1 x 10~ 4 14000 

Mn 2.24xl0~ 4 8000 

Cr 2.67xl0~ 3 16000 

Ni 7.5xl0" 3 14000 

20 P 3.1X10- 3 11000 

S 2.8xl0- 4 7500 

Gas constant R== 1.987 (cal / K • mol) 
Correction factor a: 1.2 for lC-lCr steel 

1.0 for 0.55%C steel (no cprrection) 



25 



30 

Table 3 





Symbols and physical properties in the equations of equilibrium partial pressure of CO gas 


35 


Symbol 




1C-1Cr steel 


0.55%C steel 




Pco 


Equilibrium CO gas pressure (atm) 






C 0 


Carbon content (wt%) 


1.0 


0.55 




O 0 


Oxygen content (wt%) 


0.003 


0.003 


40 


Si 0 


Si content (wt%) 


0.2 


0.2 




c L 


Carbon concentration in liquid (wt%) 








Carbon concentration in solid (wt%) 




45 


o L 

Si L 


Oxygen concentration in liquid (wt%) 
Oxygen concentration in solid (wt%) 
Si concentration in liquid (wt%) 






PS 


Density of solid ( g / cm 3 ) 


7.34 


7.30 


50 


PL 


Density of liquid ( g / cm 3 ) 


7.00 


7.00 




Kco 


Equilibrium constant ((wt%) 2 / atm) 


0.002 


0.0021 




KSiQ 2 


Equilibrium constant ((wt%) 2 / atm) 


1.94x10" 7 


7.21 xlO' 7 


55 


(Kco and are the values at the average temperatures of 1390°C and 1443°C in solidification tem- 
perature range, respectively) 




k Fe-C 


Equilibrium partition ratio in Fe-C system 


0.39 


0.37 
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Table 3 (continued) 



Symbols and physical properties in the equations of equilibrium partial pressure of CO gas 


Symbol 




IC-ICr steel 


0.55%C steel 


k Fe-Si 


Equilibrium partition ratio in Fe-O system 
Equilibrium partition ratio in Fe-Si system 
Constant in Eq. (38) 
Constant in Eq. (39) 


0.076 
0.5 


0.076 
0.5 

14.6/( Ps g s + p L g L ) 
19.5 / (ps9s + Pl9l) 


(Note: a c and a Q are obtained by applying the state equation of gas to CO gas pore ) 


ASi0 2 
Y 


Amount of Si0 2 (wt% 

Constant in Eq. (44) 


) 


0.467 



Claims 

1 . A continuous casting apparatus characterized by comprising means for exerting electromagnetic body force to 
exert electromagnetic body force (Lorentz force) toward casting direction of cast piece, in continuous casting. 

2. The continuous casting apparatus described in claim 1 , in which said means for exerting electromagnetic body 
force is characterized by exerting electromagnetic body force (Lorentz force) toward casting direction onto the solid- 
liquid coexisting zone of said cast piece. 

3. The continuous casting apparatus described in claim 1, in which said means for exerting electromagnetic body 
force is characterized by exerting electromagnetic body force (Lorentz force) toward casting direction onto the solid- 
liquid coexisting zone in the vicinity of the final solidification portion of said cast piece. 

4. The continuous casting apparatus described in claim 1, in which said means for exerting electromagnetic body 
force is characterized by exerting electromagnetic body force (Lorentz force) of the magnitude necessary to hold 
the interdendritic liquid pressure larger than the critical pressure of porosity formation onto the solid-liquid coexist- 
ing zone in the vicinity of the final solidification portion of said cast piece. 

5. The continuous casting apparatus described in claim 4, in which said means for exerting electromagnetic body 
force is characterized by having a calculating means for calculating the magnitude of said electromagnetic body 
force (Lorentz force) and the exerting region in the solid-liquid coexisting zone in the vicinity of the final solidification 
portion of said cast piece, on the basis of the operating parameters comprising at least alloy composition, cross- 
sectional shape and dimensions, casting temperature, casting speed and cooling conditions at the surface of said 
cast piece. 

6. The continuous casting apparatus described in daim 4, in which said means for exerting electromagnetic body 
force is characterized by comprising a calculating means for calculating the magnitude of said electromagnetic 
body force (Lorentz force) and the exerting region in the solid-liquid coexisting zone in the vicinity of the final solid- 
ification portion of said cast piece, on the basis of the operating parameters comprising at least alloy composition, 
cross-sectional shape and dimensions, casting temperature, casting speed, cooling conditions at the surface, dis- 
solved gas content in liquid phase and the deformation velocity due to bending, unbending, reduction of said cast 
piece. 

7. The continuous casting apparatus described either in claim 5 or 6, in which said calculating means is characterized 
by calculating the position of porosity formation, on the basis of the pressure drop of the liquid phase caused by the 
interdendritic liquid flow in the solid-liquid coexisting zone of said cast piece. 

8. The continuous casting apparatus described either in claim 5 or 6, in which said calculating means is characterized 
by comprising a correction means for correcting said magnitude of electromagnetic body force (Lorentz force) and 
said exerting region, on the basis of measured values. 

9. The continuous casting apparatus described in claim 8, in which said correction means is characterized by con- 
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ducting calculation processing, on the basis of the measured data by experiments. 

1 0. The continuous casting apparatus described in claim 8, in which said correction means is characterized by having 
a function to perform real time feed back control for said magnitude of electromagnetic body force, said exerting 

5 region and operating parameters, on the basis of the measured data of operating parameters. 

1 1 . The continuous casting apparatus described either in claim 5 or 6, in which said calculating means is characterized 
by having displaying means for displaying in real time the solidification process of said cast piece, on the basis of 
the measured values. 

10 

12. The continuous casting apparatus described in claim 1, in which said means for exerting electromagnetic body 
force is characterized by being installed concurrently with reduction means for giving reduction gradient onto said 
cast piece. 

75 13. The continuous casting apparatus described in claim 12, in which said reduction means is characterized by giving 
the reduction gradient smaller than the solidification contraction gradient, through the surface of said cast piece, in 
the solid-liquid coexisting zone in the vicinity of the final solidification portion of said cast piece. 

14. The continuous casting apparatus described in claim 12, in which said reduction means is characterized by com- 
20 prising at least a pair of rolls with said cast piece put between these rolls. 

15. The continuous casting apparatus described in claim 12, in which said reduction means is characterized by being 
composed of giving reduction force of the magnetic attractive action of said means for exerting electromagnetic 
body force. 

25 

16. The continuous casting apparatus characterized by comprising means for exerting electromagnetic body force to 
exert electromagnetic body force (Lorentz force) toward casting direction of cast piece in continuous casting, this 
means for exerting electromagnetic body force being installed with plural number of sliding electrodes attached in 
contact with the side surfaces of said cast piece, and at least one pair of superconductive coils being arranged so 

30 as to cross the direction of the current between these electrodes with the direction of generated magnetic field. 

17. The continuous casting apparatus described in claim 16, in which said means for exerting electromagnetic body 
force is characterized by comprising plural number of pairs of rolls arranged with cast piece put between them for 
supporting the cast piece or for compressing the cast piece to give larger reduction amount as it goes downstream 

35 of the cast piece. 

18. The continuous casting apparatus described in claim 17, in which said roll reduction unit is characterized by com- 
prising roll driving unit to properly balance the drawing resistant force resulting from giving reduction gradient to 
cast piece, and the electromagnetic body force exerted toward casting direction. 

40 

19. The continuous casting apparatus described either in claim 17 or 18, in which means for giving reduction force is 
characterized by using independently controllable fluid hydraulic cylinder. 

20. The continuous casting apparatus described either in claim 16 to 19, which is characterized by comprising means 
45 for removing oxidized surface layer of cast pece at the upstream side of said sliding electrodes. 

21. The continuous casting apparatus described in claim 20, which is characterized by comprising means for shielding 
the sliding part of said sliding electrodes and cast piece, and cutting part of cast piece surface with oxidation pre- 
venting gas atmosphere. 

50 

22. The continuous casting apparatus described in claim 16, in which means for exerting electromagnetic body force 
is characterized by comprising a mechanism movable forward and backward in casting direction and fixable at a 
prescribed position. 

55 23. The continuous casting apparatus characterized by comprising means for exerting electromagnetic body force to 
exert electromagnetic body forces (Lorentz forces) toward casting directions onto the cast pieces in the range of 
space which include plural number of cast pieces, in the continuous casting where these plural number of cast 
pieces are concurrently, simultaneously cast. 
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24. The continuous casting apparatus described either in claim 16 to 23, which is characterized by comprising roll 
reduction unit in the complete solid region of cast piece at a downstream side of the means for exerting electromag- 
netic body force, and thereby giving braking force created by the frictional force with the cast piece due to the roll 
reduction, onto the cast piece, so as to correspond to the electromagnetic body force exerted toward casting direc- 

5 tion. 

25. The continuous casting apparatus characterized by comprising, means for exerting electromagnetic body force to 
exert almost uniform electromagnetic body force (Lorentz force) toward the opposite direction to the casting direc- 
tion in continuous casting, said means for exerting electromagnetic body force having plural number of sliding elec- 

10 trodes being attached in contact with the both sides of said cast piece; at least one pair of superconductive coils 
being arranged so as to cross the direction of the current between the electrodes with the direction of magnetic 
field; and calculating means for calculating the porosity formation range in the solid-liquid coexisting zone in the 
vicinity of the final solidification portion of said cast piece, on the basis of the operating parameters comprising at 
least alloy composition, cross-sectional shape and dimensions, casting temperature, casting speed and cooling 

15 conditions at the surface of said cast piece. 

26. The continuous casting apparatus described in claim 25, wherein said means for exerting electromagnetic body 
force is characterized by comprising calculating means for calculating the porosity formation range in the solid-liq- 
uid coexisting zone in the vicinity of the final solidification portion of said cast piece, on the basis of the operating 

20 parameters comprising at least alloy composition, cross-sectional shape and dimensions, casting temperature, 
casting speed, cooling conditions at the surface, dissolved gas content in liquid phase and the deformation velocity 
due to bending, unbending, reduction of said cast piece. 

27. The continuous casting apparatus described either in claim 25 or 26, wherein said exerting means for electromag- 
25 netic body force is characterized by arranging plural number of pairs of rolls with cast piece put between them for 

supporting the cast piece or for compressing the cast piece to give larger reduction amount as it goes downstream 
of the cast piece. 

28. A continuous casting method characterized by comprising a step for exerting electromagnetic body force to exert 
30 electromagnetic body force (Lorentz force) toward casting direction of cast piece, in continuous casting. 

29. The continuous casting method described in claim 28, in which the exerting region of electromagnetic body force 
is characterized by being the solid-liquid coexisting zone of said cast piece. 

35 30. The continuous casting method described in claim 28, in which the exerting region of electromagnetic body force 
is characterized by being the solid-liquid coexisting zone in the vicinity of the final solidification portion of said cast 
piece. 

The continuous casting method described in claim 28, in which the exerting region of electromagnetic body force 
is characterized by being the solid-liquid coexisting zone in the vicinity of the final solidification portion of said cast 
piece, and said electromagnetic body force (Lorentz force) is characterized by being that of the magnitude neces- 
sary to hold the interdendritic liquid pressure larger than the critical pressure of porosity formation. 

The continuous casting method described in claim 31 , in which said electromagnetic body force and said exerting 
region are characterized by using the values of the magnitude and of the exerting region in the solid-liquid coexist- 
ing zone in the vicinity of the final solidification portion respectively, which are calculated on the basis of the oper- 
ating parameters comprising at least alloy composition, cross-sectional shape and dimensions, casting 
temperature, casting speed and cooling conditions at the surface of said cast piece. 

so 33. The continuous casting method described in claim 31 , in which said electromagnetic body force and said exerting 
region are characterized by using the values of the magnitude and of the exerting region in the solid-liquid coexist- 
ing zone in the vicinity of the final solidification portion respectively, which are calculated on the basis of the oper- 
ating parameters comprising at least alloy composition, cross-sectional shape and dimensions, casting 
temperature, casting speed, cooling conditions at the surface, dissolved gas content in liquid phase and the defor- 

55 mation velocity due to bending, unbending, reduction of said cast piece. 

34. The continuous casting method described either in claim 32 or 33, in which the calculation of the exerting zone of 
said electromagnetic body force (Lorentz force) is characterized by being obtained from the position of porosity for- 
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mation which is obtained on the basis of the pressure drop of the liquid phase caused by the interdendritic liquid 
flow in the solid-liquid coexisting zone of said cast piece . 

35. The continuous casting method described either in claim 32 or 33, in which the magnitude of said electromagnetic 
5 body force (Lorentz force) and said exerting region are characterized by being corrected by the corrected values 

obtained on the basis of the measured values. 

36. The continuous casting method described in claim 35, in which the said corrected values are characterized by 
being obtained on the basis of the measured data by experiments. 

10 

37. The continuous casting method described in claim 28, in which reduction gradient is characterized by being given 
in said exerting region of electromagnetic body force or its vicinity. 

38. The continuous casting method described in claim 37, in which said reduction gradient is characterized by being 
15 smaller than the solidification contraction gradient in the vicinity of the final solidification portion of said cast piece, 

and being given through the surface of said cast piece. 
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Figure 2 
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Figure 3 (a) Equilibrium phase diagram 




(b) Solute distribution for an equilibrium solidification 
type alloy at the temperature T 
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(C) Solute distribution for a nonequil ibrium solidification 
type alloy at the temperature T 
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Figure 10(a) Main program 
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Figure 10(b) 

Analysis of the flow field by the extended SIMPLER method 
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Figure 11 ( a ) 
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Figure 12 
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Figure 14 



1500 



1400 - 



T(C) 



1300 - 



1200 



1100 




Figure 15 



(a) lC-lCr Bearing steel 



1453 C 



T(C) 




0.2 0.4 0.6 0.8 

Vol. fraction of solid gs 



l.o 




83 

BNSDOCtD: <EP 090959BA1_J_> 



EP 0 909 598 A1 



Figure 16 
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Figure 18 
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Figure 20 
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Figure 21 
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Figure 22 
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Figure 23 

(a) Schematic diagram (b) Horizontal cross section 
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Figure 24 
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Figure 26 
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Figure 27 
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Figure 28 
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Figure 29 
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Figure 30 
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Figure 31 
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Figure 32 
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Figure 34 
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Figure 35 
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Figure 36 
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Figure 37 
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Figure 39 
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Figure 41 
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Figure 42 
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Figure 43 
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Figure 44 
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Figure 45 
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Figure 46 
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Figure 47 
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Figure 49 
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Figure 50 
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Figure 52 
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Figure 53 
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Figure 55 
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Figure 57 
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Figure 59 

(a) BB cross section of Figure 58(a) 
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Figure 62 
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Figure 64 
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